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ABSTRACT 


The  changes  in  VLF  radio  noise  accompanying  a  solar  flare  and  its  asso¬ 
ciated  sudden  ionospheric  disturbance  (SID)  are  surveyed.  The  onset  time  of  the 
effects  is  about  10  minutes,  but  may  be  less  than  2  minutes  for  perhaps  5  per¬ 
cent  of  the  incidents.  The  alteration  in  propagation  conditions  implies  an  in¬ 
crease  in  noise  on  frequencies  exceeding  20  kc  and  a  decrease  below  10  kc.  If 
a  local  thunderstorm  develops,  the  time  phasing  is  very  similar  to  that  of  an 
SID  incident,  but  noise  at  all  frequencies  is  enhanced. 

A  detailed  examination  of  the  formative  processes  of  the  lower  ionosphere 
is  made.  Values  of  the  rate  coefficients  for  attachment,  detachment  (photo, 
collisional,  and  associative),  electron-ion  recombination,  and  ion-ion  recom¬ 
bination,  are  deduced  for  heights  from  40  to  100  km.  Over  the  same  height 
range,  the  ionizing  contributions  due  to  the  major  influences  of  cosmic  rays, 
solar  X  rays,  and  solar  Lyman  o?  and  to  certain  subsidiary  causes  arc  examined. 
The  balance  equations  for  the  lower  ionosphere  are  solved  for  the  cases  of 
equilibrium  conditions  and  of  the  departure  from  equilibrium  associated  with 
the  approach  of  night.  Two  specific  examples  indicate  that  contours  of  electron 
density,  and  therefore  of  ionospheric  conductivity,  are  far  from  parallel  to  the 
surface  of  the  earth.  It  follows  that  the  attenuation  coefficient  for  VLF  propa¬ 
gation  by  day  can  differ  according  to  the  solar  zenith  angle  by  as  much  as  20 
percent  from  the  mean  daytime  value. 

A  model  of  the  solar  X-ray  flux  during  a  flare  is  derived  and  its  ionizing 
effects  examined.  The  flux  peaks  in  5  minutes  but  it  is  found  that  the  maximum 
ionospheric  effect  lags  behind,  especially  at  the  lower  heiglits.  Suii)risingly 
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enough,  the  persistence  of  the  added  ionization  does  not  increase  moiiotonieally 
with  increasing  height. 

The  ionization  produced  by  an  X-ray  burst  from  a  nuclear  explosion  in 
space  is  deduced.  Differences  in  Ibe  onset  phase,  the  spectral  distribution  of 
the  X  rays,  and  the  zenith-angle  dependence,  as  compared  with  a  natural  SID, 
are  demonstrated.  The  decay  of  the  additional  ionization  is  exam.ined  and  it  is 
shown  that  this  often  involves  three  phases;  the  establishment  of  an  attachment- 
detachment  balance,  a  plateau  with  little  decay,  and  a  final  decrease  by  recom¬ 
bination.  An  example  indicates  that  if  the  same  amount  of  extra  ionization  is 
introduced  at  heights  of  70,  75,  and  80  km,  then,  after  about  an  hour,  more  of 
the  additional  ionization  remains  at  70  km  than  at  either  of  the  other  heights. 
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1.  INTRODUC  riON 


The  very-low-frequency  (VLF)  band  in  the  radio  spectrum  is  officially 
defined  as  covering  the  frequency  range  of  3  to  30  kc.  Radio  waves  of  these 
frequencies  are  normally  propagated  in  the  quasi -waveguide  represented  by  the 
earth  and  the  lower  ionosphere.  The  special  case  of  the  whistler  mode  [Hclliwell 
and  Morgan  (1959)],*  in  which  the  VLF  signals  are  guided  along  a  magnetic  line 
of  force  and  penetrate  the  ionosphere,  will  not  be  considered  in  this  report. 

Experimental  work  upon  VLF  propagation  is  reasonably  consistent  in 
suggesting  that  sky  waves  are  reflected  from  a  height  of  about  70  km  by  day; 
at  night  the  reflection  height  is  some  90  km  |  Braeewell  et  al.  (1951)].  Also,  the 
theoretical  approach  in  which  the  ionosphere  is  considered  to  be  a  homogeneous 
conductor  sharply  bounded  at  the  heights  quoted  above  has  been  effective  in 
explaining  many  of  the  major  features  experimentally  observed  in  VLF  propagation 
[Wait  (1957,  I960)],  The  model  is  that  of  a  waveguide  in  which  the  walls  — 
spherically  concentric — are  formed  by  the  earth  and  the  lower  ionosphere;  both 
of  the  walls  are  of  uniform  finite  conductivity.  Later  theoretical  work  by  Wait 
(1958)  suggests  that  a  better  agreement  with  experimental  results  is  found  if  an 
increase  of  ionization,  and  hence  of  conductivity,  with  penetration  into  the  iono¬ 
spheric  boundary  is  postulated.  However,  it  is  still  true  that  the  dominant  factor 
in  VLF  propagation,  especially  by  day,  is  the  ionization  at  heights  below  100  km. 

The  lower  ionosphere,  which  may  be  roughly  identified  with  the  D  region,  is 
complicated  and  far  from  being  entirely  understood.  The  principal  iomzing 
agencies  are  cosmic  rays  and  the  solar  X-ray  and  Lyman-O'  (wavelength  1216A) 
radiations;  the  latter  arc,  of  course,  ab.scnt  at  night.  Subsidiary  ionization 


‘References  arc  given  at  the  end  of  each  section. 
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results  from  the  action  of  meteors,  of  other  solar  ultraviolet  radiations,  and  of 
incoming  particles  with  origins  in  the  sun  or  in  the  van  Allen  belts  around  the 
earth.  The  immediate  action  of  the  ionizing  influences  is  to  produce  free 
electrons  and  positive  ions  of  molecular  or--above  about  90  km — atomic  size. 
Subsequently,  the  electrons  may  become  attached  to  neutral  particles  thus 
forming  negative  ions;  the  reverse  process  of  detachment  may  operate:  and  there 
will  be  recombination  between  electrons  and  positive  ions,  and  also  between 
positive  and  negative  ions.  Under  quiet  solar  conditions,  balance  equations,  in 
which  all  the  above  processes  are  represented,  can  be  formulated  and  solved. 

There  are,  however,  many  uncertainties  in  the  coefficients  involved.  Also,  it 
is  quite  possible  that  ions  of  larger  size  than  single  molecules  are  present  in  the 
lower  ionosphere,  a  fact  that  would  further  complicate  the  balance  equations.  It 
is  known  that  above  90  km  the  ions  arc  of  molecular  or  atomic  character;  at  30 
km  the  dominant  ion  is  the  small  or  cluster  ion  of  atmospheric  electricity  con¬ 
sisting  of  a  single  molecular  ion  surrounded  by  a  group  of  perhaps  4  to  12 
neutral  molecules;  the  height  at  which  the  transition  from  cluster  ions  to  molecular 
ions  takes  plaec  remains  uncertain. 

When  a  flare  is  observed  upon  the  sun,  normally  in  the  light  of  the  hydrogen 
lla  line  (6563A),  a  sudden  ionospheric  disturbance  (SID)  usually  follows  within  a 
few  minutes.  The  SID  is  almost  certainly  due  to  an  enhancement,  associated 
with  the  flare,  especially  of  the  solar  X-ray  flux  [  Friedman  (1959)).  Additional 
ionization  is  produced  in  the  ionosphere,  this  being  detectable  in  several  ways: 
for  e.xample,  a  sharp  kick  known  as  a  crochet  occurs  upon  records  of  the  earth'  s 
magnetic  field.  There  are  also  pronounced  absorption  effects  upon  radiowaves 
passing  through  the  D  region;  this  is  because  the  absorption  coefficient  for  such 
nondeviating  waves  is  proportional  to  the  product  of  electron  density,  N,  and 
electron  collisional  frequency,  r,  so  that  the  increase  in  N  represents  an 
intensified  absorption.  The  absorption  is  usually  detected  in  two  ways.  Sources 
of  cosmic  noise  may  be  monitored,  normally  at  a  frequency  (about  20  Me) 
slightly  above  the  lower  limit  that  would  penetrate  the  F  layer:  the  drop  in 


signal  accompanying  an  SID  is  then  often  rcfcri'cd  to  as  a  sudden  cosmic  noise 
absorption  (SCA  or  SCNA).  Alternatively,  signals  generated  by  ground-based 
transmitters,  at  frequencies  (3-15  Me)  that  arc  usually  reflected  from  one  of 
the  ionospheric  layers  above  D  (E  or  F),  may  be  observed.  In  this  case  there 
is  a  multiple  path  through  the  D  region  and  the  absorption  produces  a  short¬ 
wave  fade-out  (SWF), 

The  effects  of  an  SID  upon  VLF  propagation  arc  also  marked.  Here  the 
influence  of  the  inereased  ionization  is  to  depress  the  lower  boundary  of  the 
D  region  and  thus  to  narrow  the  earth-ionosphere  waveguide  from  its  normal 
daytime  dimension  of  70  km.  At  the  same  time,  the  conductivity  of  the  iono¬ 
sphere  may  be  enhaneed.  There  are  two  principal  SID  phenomena  detectable 
at  VLF.  The  first  of  these  is  the  sudden  phase  anomaly  (SPA);  when  the  phase 
of  a  long-wave  transmitter  is  continuously  monitored,  a  rapid  change,  resulting 
Irom  the  lowering  of  the  reflcetion  points  for  sky  waves,  is  observed.  Secondly, 
an  SID  alters  the  propagation  characteristic  of  the  earth-ionosphere  waveguide 
so  that  the  attenuation  of  some  frequencies  in  the  VLF  band  (3-30  kc)  is  increased 
while  that  of  others  is  diminished.  The  changes  arc  most  conveniently  detected 
l)y  recording  the  radio-noise  or  VLF  atmospherics  generated  by  lightning:  this 
source  radiates  at  all  frequencies  in  the  VLF  band  so  that  the  influence  of  the 
SID  can  be  determined  us  a  function  of  frequency.  The  phenomenon  is  often 
referred  to  as  a  sudden  enhancement  of  atmospherics  (SEA),  although  this  de¬ 
scription  is  somewhat  misleading:  the  enhancement  is  only  at  certain  frequencies 
while  at  others  there  is  a  decrease  in  signal  strength. 

This  report  rc()rcscnts  Part  I  of  the  Final  Report  under  U.  S.  Air  Force 
Contract  AF  49(638)-1081.  The  complete  Final  Report  will  be  concerned  with 
various  aspects  of  the  interrelation  between  an  SID,  cither  of  natural  or  man¬ 
made  origin,  and  the  associated  changes  in  attenuation  coefficients  for  propaga¬ 
tion  at  frequencies  within  the  VLF  band.  The  SEA  phenomenon  will  receive 
jirincipal  attention.  This  report  (Part  1)  considers  two  main  subjects.  Fii'st, 
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the  phenomenology  of  a  natural  SID  and  its  associated  changes  in  VLF  radio 
noise  are  discussed.  Second,  a  model  of  the  lower  ionosphere  and  of  its  changes 
under  various  conditions  is  derived;  in  particular,  variations  during  an  SID  and 
following  a  nuclear  explosion  in  space  are  considered.  The  application  to  VLF 
propagation  is  examined.  Part  II  of  the  Final  Report  will  mainly  deal  with  the 
experimental  data  obtained  from  an  analysis  of  whistler  tapes,  and  with  the 
development  of  simple  instruments  for  monitoring  the  changes  in  VLF  radio 
noise  accompanying  an  SID. 
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2.  PHENOMENOLOGY  OF  SUDDEN  lONOSPllElUC  DISTURBANCES 
AND  THE  ASSOCIATED  CHANGES  IN  YLE  RADIO  NOISE 


2.1  INTRODUCTION 

The  solar  ()ull)urst,  known  as  a  Hare,  usually  occurs  in  association  wilh 
an  active  spot-f^roup  upon  tlic  sun.  The  tlari'  rep:ion,  as  observed  in  Ho  li<;ht, 
brifrbtens  and  increases  in  area  for  a  few  minutes;  sul>sec|iiontly  there'  is  a 
return  to  nornml  which  may  occupy  perhaps  half  an  hour.  Flares  art'  classifit'd 
primarily  in  terms  of  area,  and  the  basic  catef;oi’ics  in  oriltM’  of  inci'casin^' 
importance  art'  flenotetl  by  1.  2,  anti  The  supitlementary  tlesifrnatit)ns  1-,  1  i. 

2+,  anti  :d  are  witlely  employetl.  AllhouRh  tht'  classification  is  basetl  uptm 
area,  increasing  flare  impt>rtancc  also  implies  an  increase  in  total  tluratittn. 
in  the  line-witith  of  lit)  at  the  rna.ximum,  anti  in  the  intensity  ttf  llo  iTom  thi' 
flare  rcgit)n  as  comparctl  with  the  continuum  frtim  atljacent  art'as.  Tht'st' 
features  are  illustratcfl  in  Fig.  2.1.  When  flart'S  are  divitlt'il  into  Categtirit's 
1,  2,  anti  :i,  it  is  ft)untl  that  72  percent  of  the  incitlcnts  art'  ttf  tlie  first  imporlanct' ; 
2r)  iK'rcent  are  in  Category  2:  anti  only  2,  percent  are  elassiiit'd  as  largt'  Hart's. 

Hti  light  penetrates  to  the  surfact'  t)f  the  earth  withttul  effeclivt'  ahstirption, 
anti  thus  enables  the  tlevelopment  of  a  flare  to  be  rcatlily  monitoretl.  llowt'vt'r. 
any  stHar  radiations  that  influence  the  ionttsithcre  must  suffer  apprt'ciablt' 
absorption  there:  it  follows  that  such  ratliations  are  not  t'asily  rlf'tt'ctahlt' , 
except  in  the  high  atmosphere.  Divt'ct  informatitm  uptm  their  variation  tluring 
an  SID  must  ix’  largely  obtainetl  frttm  rocket  antl--t'specially--satellitt' 
obser\ations.  Little  such  information  is  iiresently  available  Inft'rt'nces  frttm 
itmosphcric  phenomena  still  supply  most  t)f  the  itnowh'tlgt'  regartling  Hit'  behavior 
tif  solar  ionizing  ratliations  tim  ing  flares. 
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IMPORTANCE  OF  FLARE 


FIG.  2.1  SOME  FEATURES  OF  SOLAR  FLARES 


For  ionization  to  be  produced  somevvliere  in  llie  ionosphere,  the  ionizing 
radiation  must  liave  a  photon  energy  >9.5  ev,  tJie  energy  representing  ionization 
of  the  most  easily  ionizable  moleeule,  NO.  In  terms  of  wavelengtli,  X  ,  this  im¬ 
plies  a  wavelengtli  of  less  than  some  1300 A.  Within  this  range  the  harder  X  rays 
{A<8a)  and  Lyman  -o  radiation  (X=  1216A)  are  capable  of  producing  ionization 
well  below  100  km;  the  intermediate  X-ray  and  ultraviolet  radiations  are  mostly 
absorbed  in  the  E  and  F  regions.  Direct  evidence  from  the  NHL  rocket  program 
[  Chubb,  Friedman  and  Kreplin,  {I960)  1  has  shown  that  there  is  certainly  an  en¬ 
hancement  of  the  solar  X-ray  flux  during  a  flare.  The  increase  extends  over  a 
range  of  X-ray  wavelengths  capable  of  giving  ionization  in  both  the  D  and  E 
regions.  Also  the  NRL  rocket  program  has  failed  to  indicate  any  appreciable 
change,  associated  with  a  flare,  in  the  lyman  -o  flux.  No  direct  evidence  is 
available  for  the  wavelength  range  between  Lyman  n  and  X  ray,  but  tlie  obser¬ 
vation,  at  HF,  of  flare-induced  frequency  shifts,  by  Kanellakos,  Chan,  and 
Villard  (1962),  strongly  suggests  that  e.xtra  ionization  is  produced  at  heights 
between  120  and  200  km.  This  in  turn  would  imply  an  increase  in  the  ultraviolet 
flux  in  the  wavelength  range  of  100  to  1000  A. 

There  is  little  evidence  for  the  ph  .sing  of  the  ionizing  radiations  with  regard 
to  Ho  during  a  flare.  One  of  the  few  direct  observations  was  obtained  by  the  solar 
satellite  SR  1  [Kreplin,  Chubb,  a.nd  Friedman  (1962)],  for  a  flare  on  August  6, 
1960.  This  indicated  that  the  X-ray  (<8A)  flux  started  to  rise  at  about  1510  U.T. 
reaching  a  maximum  at  approximately  1513;  the  Hn  start  was  at  1506  and  the 
maximum  at  1514;  and  the  Lyman  -o  fliLX  remained  constant  between  1509  and 
1514.  The  work  of  Kanellakos  et  al.  (1962)  shows  that  the  maximum  HF  shift 
occurs  1  to  4  minutes  before  the  maximum  phase  in  Ho'.  Generally,  the  avail¬ 
able  observations  and  inferences  are  consistent  in  suggesting  that  any  ultraviolcl 
radiation  is  likely  to  precede  the  Ho,  but  that  tlie  X-ray  and  Ho  phasing  almost 
coincide.  The  hardest  (highest-frequeney)  X  rays  are  most  pronounced  when  the 
total  X-ray  intensity  is  maximum;  a  similar  hardening  to  maximum  with  a  sub¬ 
sequent  softening  jirobably  also  holds  for  the  ultraviolet  flux. 


The  ionospheric  consequences  of  the  ionizing  radiations  depend,  of  course, 
upon  their  relative  penetration.  The  VLF  effects,  namely  the  SPA  and  SEA,  are 
controlled  by  hard  X  radiation  and  the  consequent  increase  in  ionization  below 
about  70  km.  Phase  and  frequency  changes  at  megacycle  frequencies,  on  the 
other  hand,  arc  caused  by  ionization  at  considerably  higher  levels,  perhaps 
even  above  150  km  [Morriss  (I960)];  this  is  presumably  due  to  ultraviolet 
radiation.  Absorption  is  an  integrated  effect  throughout  much  of  the  ionosphere 
and  in  consequence  the  SWF  and  SCNA  phenomena  cannot  be  easily  related  to 
increases  in  ionization  at  specific  heights. 

Since  a  broad  range  of  ionizing  radiations,  and  consequently  of  ionospheric 
phenomena,  is  involved  in  a  solar  flare,  the  various  effects  are,  inevitably,  not 
very  closely  related.  Although  there  is  an  interconnection  in  the  average 
behavior,  as  shown,  for  example,  in  the  straight-line  relationship  postulated 
by  Shain  and  Mitra  (1953)  between  SPA  and  SCNA  magnitudes,  there  are  wide 
differences  between  the  individual  cases.  Some  flares,  for  instance,  give 
strong  VLF  effects,  but  minor  absorption  changes;  others  may  show  considerable 
absorption,  but  little  alteration  in  VLF  propagation.  This  variability  cannot  be 
overstressed;  it  is  entirely  incorrect  to  assume  that  the  relative  importance  of 
the  ionospheric  effects  associated  with  a  solar  flare  is  rigidly  fixed. 

2.2  SOME  SID  STATISTICS 

The  spread  in  the  phenomena  accompanying  a  flare  has  been  emphasized 
in  the  preceding  section;  nevertheless,  the  average  behavior  can  still  be 
instructive,  particularly  if  the  statistical  variability  incorporated  in  the 
averaging  procedure  is  also  considered.  Accordingly,  the  six-month  period 
of  January -June  1959  was  selected  and  the  information  available  in  "Compilations 
of  Solar-Geophysical  Data,"  CRPL-F,  Part  B,  issued  by  the  National  Bureau 
of  Standards,  was  examined.  Over  a  hundred  incidents  were  chosen  for  which 
SCNA  and  SEA  statistics  (at  27  kc)  were  available,  the  SEA  being  regarded  as 


representative  of  VLF  effects  and  the  SCNA  bein;;'  considered  as  indicative  of 
absorption;  in  about  50  percent  of  cases,  llfi  flare  information  could  also  be 
definitely  associated  with  the  ionospheric  incident.  The  statistics  examined 
eonsisted  of  duration  for  SCNA,  SEA,  and  Hcj;  flare;  the  rise  time  (time  from 
start  to  maximum)  for  the  same  three  phenomena;  the  delay  in  the  start  of  the 
SEA  and  SCNA  as  compared  with  that  for  Ha;  and  the  delay  in  the  maximum  of 
the  SEA  and  SCNA  again  as  compared  with  that  for  Ha. 

The  results  obtained  are  shown  in  Figs.  2.2  to  2.6.  They  are  in  fair 
agreement  (see  Table  2.1)  with  some  of  the  previous  work  of  Ellison  (1957, 
1960)  who  tabulated  SWF  and  SEA  statistics.  The  agreement  is  best  for  the 
two  sets  of  SEA  data,  and,  in  view  of  the  fact  that  both  the  SWF  and  SCNA  are 
basically  absorption  effects,  some  of  the  differenees  indicated  between 
Ellison's  results  for  SWFs  and  the  present  results  for  SCNAs  are  suprising. 

Certain  features  of  the  data  are  worth  mentioning.  It  seems  likely  that 
the  peak  in  ionospherie  phenomena  will  lag  somewhat  behind  the  maximum  of 
the  ioni/'.ing  radiation,  sinee  until  equilibrium  is  established  between  produetion 
and  dissipation  of  eleetrons,  the  effects  will  be  cumulative.  It  is  certain  th.at 
upon  eessation  of  the  ionizing  flux  the  ionosphere  will  take  an  appreciable 
time--depending  upon  the  height  eoncerned--to  return  to  normal.  Many  of  the 
SBD  incidents  examined  show  the  Ha  flare  commeneing  before  either  the  SEA 
or  the  SCNA,  the  maximum  in  the  Ha  flux  preceding  that  of  the  SEA  and  SCNA, 
and  the  Ha  light  returning  to  normal  well  before  complete  recovery  of  the 
ionosphere.  Such  incidents  are  entirely  compatible  with  the  concept  that  the 
variation  of  the  ionizing  radiation  is  closely  related  to  that  of  Ha.  However, 
Fig.  2.4  indicates  that  in  a  small  but  significant  proportion  of  instances  the 
onset  of  an  SCNA  or  SEA  precedes  that  of  the  Ha  flare.  Furthermore,  it  is  not 
uncommon  for  the  SEA  and  SCNA  effects  to  have  ceased  while  the  Ha  emission 
is  still  active;  in  some  instances  the  duration  of  the  Ha  flare  may  exceed  that 
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FIG.  2.2  DURATION  STATISTICS  FOR  SID  INCIDENTS 
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FIG.  2.3  RISE  TIME  STATISTICS  FOR  SID  INCIDENTS 
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FIG.  2.4  STARTING  DELAY  STATISTICS  FOR  SID  INCIDENTS 
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FIG.  2.5  STATISTICS  FOR  DELAY  OF  MAXIMUM  FOR  SID  INCIDENTS 
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FIG.  2.6  AVERAGE  PHASING  OF  Ha  FLARE,  SCNA,  AND  SEA,  DURING  AN  SID  INCIDENT 

of  the  SEA  and  SCNA  by  60  minutes  or  more.  All  these  considerations  confirm 
the  conelusion  that  although  the  Ho-  behavior  is  a  guide  to  that  of  the  ionizing 
radiations,  there  is  no  exact  correspondence. 

One  aspect  of  SID  phenomena  that  could  bo  instructive  as  regards  tlie 
ionizing  radiation  involved  is  the  dependence  upon  solar  zenith  angle.  Accordingly, 
the  SCNA  incidents  for  the  period  March  1958  to  September  1961  were  examined, 
the  data  being  obtained  as  before  from  the  solar-geophysical  compilations  of  the 
NBS.  For  each  occurrence,  the  solar  zenith  angle  at  the  time  of  the  SCNA  and 
the  percentage  of  absorption  were  obtained  for  the  reporting  station,  and  a  plot 
(Fig.  2.7)  made  of  the  information.  (One  point  was  obtained  for  each  SCNA.)  It 
was  anticipated  that  a  trend  toward  increasing  absorption  with  decreasing  zenith 
angle  would  be  detected;  no  such  trend  is  apparent  on  Fig.  2.7,  cither  in  the 
individual  points  or  in  the  lO-degree  averages.  The  reason  may  be  that  present 
SCNA  techniques  involve  threshold  levels  that  tend  to  mask  true  variations. 
Perhaps  the  best  method  of  ascertaining  the  zenith  angle  dependence  in  a  SCNA 
would  be  to  employ  results  for  the  same  SCNA  from  several  widely  dispersed 
stations. 
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RESULTS  FOR  MARCH  1956  -  SEPTEMBER  I96i 

POINTS  AVERAGED  IN  TEN  DEGREE  INTERVALS 


2.3  THE  BEHAVIOR  DURING  AN  SEA 


2.3.1  Introduction 

Almost  all  the  routine  information  upon  an  SEA  has  been  derived 
from  recordings  at  a  frequency  of  27  kc.  Observatories  in  most  countries 
employ  an  instrument  tuned  to  this  frequency  and  integrating  the  received  atmos¬ 
pheric  noise  over  a  period  of  a  few  seconds.  The  enhancements  arc  reported 
upon  a  regular  basis  and  the  results,  which  involve  world-wide  coverage,  arc 
reproduced  in  publications  such  as  "Compilations  of  Solar -Geophysical  Data" 
of  NBS.  Information  upon  the  behavior  of  SEAs  as  a  function  of  frequency  is 
scanty.  Some  data  arc  available  from  research  investigations,  and  the  results 
from  Japan,  where  it  is  common  practice  to  record  atmospheric  noise  continu¬ 
ously  upon  several  frequencies,  arc  especially  valuable. 

2.3.2  The  Rise  Time  of  an  SEA 

The  rise  times,  that  is,  the  time  from  onset  to  maximum,  of  the 
SE.As  tabulated  in  the  NBS  publications  have  been  examined.  In  these  tabulations 
each  SEA  is  classed  in  terms  of  magnitude  according  to  the  1,  2,  3,  scale  with 
intermediate  +  or  -  signs.  Ignoring  the  latter.  Fig.  2.8  shows  histograms  giving 
the  distribution  of  the  results  according  to  class.  The  information  covers  the 
period  from  March  1958  to  February  1962  and  includes  896  SEA  incidents.  The 
mean  value  of  the  rise  time  increases  with  class :  for  SEA  incidents  of  importance 
1,  2,  and  3,  the  respective  average  rise  times  are  10.0,  12.5,  and  14.8  minutes. 
The  distributions  arc  not  symmetrical  but  arc  skewed  in  the  sense  that  small 
values  less  than  the  mean  tend  to  be  more  common  than  large  rise  times.  The 
asymmetry  seems  to  be  most  pronounced  for  the  smaller  incidents,  but  this  may 
merely  reflect  the  limited  number  of  large  SEAs.  A  normal  distribution 
centered  upon  the  mean  has  been  plotted  upon  each  of  the  histograms  ;  this  assists 
in  emphasizing  the  skewness  of  the  actual  observations.  The  standard  deviations 
appropriate  to  Classes  1,  2,  and  3  arc  5.5,  6.8,  and  7.6  minutes,  respectively. 
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Poisson  distributions  are  also  given  on  Fig.  2,8;  these  might  be  expected  to  fit 
the  data  better  than  normal  curves,  but  the  agreement  is  not  good.  For  all  the 
89G  incidents  examined,  only  11  have  rise  times  of  a  minute  or  less  while  there 
are  19  SEAs  with  rise  times  between  1  and  2  minutes.  Thus  the  chance  of  an 
SEA  having  a  rise  time  equal  to  a  minute  or  smaller  is  1.2  percent:  for  less 
than  or  equal  to  2  minutes  the  chance  is  3.4  percent.  These  statistics  are  for 
a  4-year  period  of  high  solar  activity.  It  may  thus  be  confidently  stated  that 
taken  as  a  whole  over  the  complete  solar  cycle  of  11  years,  SEA  incidents  with 
rise  times  of  less  than  2  minutes  are  unlikely  to  occur  more  often  than  seven 
per  annum. 

Kamada  (1960)  has  examined  SEA  data,  obtained  in  Japan  over  the 
5  years  of  1955  to  1959,  on  the  three  frequsneies  of  10,  21,  and  27  kc.  For  each 
frequency  the  distribution  of  the  rise  times  is  skewed  towards  the  smaller  values 
in  a  manner  somewhat  similar  to  that  of  Fig.  2.8.  However,  the  average  rise 
times  are  somewhat  less  than  those  given  in  the  preceding  paragraph.  For  the 
frequencies  of  10,  21,  and  27  kc,  the  mean  rise  times  are,  respectively,  7.8, 

7.9  and  7.4  minutes.  Also,  the  proportions  of  small  rise  times  are  rather  more 
than  those  previously  quoted.  For  rise  times  less  or  equal  to  1  minute,  and  less 
or  equal  to  2  minutes,  the  respective  percentages  are  2.3  and  13.0. 

2.3.3  Behavior  as  a  Function  of  Frequency 

The  vertical  electric  field,  E  (mv/m),  at  a  great-circle  distance  d 
(in  km)  from  a  source  radiating  P  kw  is  given  approximately  [Wait  (1958)]  by 

j  n  i  ’ 

h  I  n  .’i  in  il  n  \ 

where  a  is  the  radius  of  the  earth  (km),  h  the 
layer  (km),  \  the  wavelength  concerned  (km). 


2 

rxp{~acl)  (2.1) 

height  of  the  ionospheric  reflecting 
and  O’  the  attenuation  factor. 
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Equation  (2.1)  asRumos  the  ilominancc  of  a  sinfrle  mode  which  is  appi'oxi matcly 
true  for  (I  >  1000  km  and  frequencies  in  the  VLF  range.  The  application  to  VLF 
atmospheric  noise  is  therefore  reasonably  valid.  It  follows  that  if  the  suffix  N 
represents  normal  conditions  and  S  the  SID,  then 


(a^.  -  a.,)  2  X  10  *  lofr 

(I 


I 

</ 


10'  /o;;  (l.OH  K  ^  T  ^). 


(2.2) 


In  Eq,  (2.2)  the  units  for  o  ^  and  o  ^  are  db/1000  km,  and  the  fields  E^,  and  E 
are  taken  as  being  eompared  at  a  definite  frec|ucncy.  The  factor  h.^, /h^  =  1.08; 
this  is  the  average  value  given  by  Braeewcll  and  Straker  (1940). 

Pierce-  (19(>1)  has  applied  Eq.  (2.2)  to  various  sots  of  research  data 
in  order  to  obtain  the  variation  of  (a-,,  -  with  freciucncy.  The  three  sources 
of  information  originally  considered  by  Pierce  were  Gardner  (1950),  Obayashi 
(1960),  and  Obayashi,  Fujii,  and  Kidokoro  (1959).  The  corresponding  results 
arc  shown  on  Fig.  2.9  together  with  an  additional  curve  deduced  from  the  work 
of  Bowe  (1951).  It  is  evident  from  Fig.  2.9  that  at  frequencies  of  20  kc  and 
greater  (oq.  _  is  negative;  in  other  words,  the  attenuation  is  less  during 
an  SID  than  under  normal  conditions,  and  there  is  an  enhancement  of  atmospheric 
noise.  For  frequencies  of  10  kc  and  less,  (nq.  -  o  .^, )  is  positive,  and  atmospheric 
noise  will  therefore  decrease  in  intensity  with  the  onset  of  SID  conditions.  The 
range  from  10  to  20  kc  represents  the  intermediate  zone  where  the  transition 
from  a  decrease  to  an  enhancement  occurs. 
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FIG.  2.9  DIFFERENCE  BETWEEN  VLF  ATTENUATION  COEFFICIENTS  FOR  NORMAL 
DAYTIME  AND  FOR  SID  CONDITIONS 


The  extensive  observations  of  Kamada  (19G0,  1961)  on  the  three 
frequencies  of  10,  21,  and  27  kc  show  that  when  a  weli-authcnticatcd  SID  occurs 
there  is  an  increase  of  noise  on  21  and  27  kc  and  a  decrease  at  10  kc.  This 
result  is  entirely  consistent  with  Fig.  2.9.  Kamada  estimates  that  the  transitional 
or  turnover  frequency,  which  is  unaffected  fluring  a  SEA,  lies  Ijetween  14  and 
18  kc.  This  suggestion  is  again  in  agreement  with  Fig.  2.9,  liut  it  should  Ijc 
emphasi?,ed  that  Kamada's  estimate  is  Ijascd  on  interpolation,  and  also  that  the 
amount  of  data  represented  by  Fig,  2,9  is  very  limited.  There  are  some  indica¬ 
tions  that  the  transitional  frequency  may  vary  between  individual  flares  in  a 
manner  apparently  unrelated  to  other  characicristics  or  to  associated  ionospheric 
effects.  More  data  arc  necessary  before  such  cpiestions  can  be  entirely  clai  ifiofl. 


20 


2.3.4  Local  Thuiidorslorms 


An  SID  produces  an  SEA  on  a  frequency  exceeding  20  kc  by  changing 
propagation  conditions,  hence  reducing  the  value  of  a  in  Eq.  (2.1),  and  conse¬ 
quently  increasing  E.  Obviously  alterations  in  the  other  factors  in  Eq.  (2.1)  can 
also  have  the  effect  of  enhancing  E;  such  alterations  are  due  to  changes  in  the 
level  of  thunderstorm  activity.  For  instance,  if  existing  storms  increase  in  in¬ 
tensity,  this  implies  a  greater  value  of  P  and  hence  an  enhancement  in  E;  if  the 
storms  move  closer  to  the  recording  station,  this  gives  a  decrease  in  the  effec¬ 
tive  d  and  again  an  increase  in  E;  if  new  thunderstorms  develop,  particularly  at 
short  ranges,  once  more  E  is  enhanced.  In  the  absence  of  other  information, 
observations  of  atmospheric  noise,  at  the  frequency  of  27  kc  usually  employed, 
cannot  distinguish  between  an  SEA  associated  with  an  SID  and  one  due  to  changes 
in  thundery  activity,  unless  there  are  characteristic  differences  between  the  two 
types  of  SEA. 

The  frequency  dependence  enables  the  two  types  of  SEA  to  be  readily 
distinguished.  Section  2.3.3  indicates  that,  accompanying  an  SID,  a  reciprocal 
effect  applies;  there  is  an  enhancement  of  noise  at  some  frequencies  and  a  de¬ 
crease  at  others.  Changes  in  the  level  of  thundery  activity  will,  however,  pro¬ 
duce  uniform  behavior;  at  all  frequencies,  noise  will  either  be  increased  or 
diminished. 

Japanese  observations  for  the  year  1958  [  IGY  Data  (1959)]  are  avail¬ 
able,  the  sudden  changes  in  atmospheric  noise  observed  at  frequencies  of  10,  21, 
and  27  kc  being  listed.  By  means  of  the  criterion  of  reciprocal  or  similar  ampli¬ 
tude  behavior  at  10  and  27  kc,  it  is  easy  to  separate  the  SEA  incidents  associated 
with  solar  disturbances  from  those  probably  due  to  changed  thunderstorm  activity. 
Various  statistics  derived  for  the  two  types  of  SEA  are  shown  on  Table  2.2. 
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TABLE  2.2 


COMPARISON  OF  AVERAGE  CHARACTERISTICS  FOR  SEAS  ASSOCIATED 
RESPECTIVELY  WITH  SIDS  AND  CHANGES  OF  THUNDERY  ACTIVITY 


Mean  Time*  That 

Mean  Rise  Time* 

Mean  Duration* 

Onset  at  10  kc  Precedes 

'rypc  of  SEA 

10  kc 

27  kc 

10  kc 

27  kc 

That  at  27  kc 

Associated 
with  SID 

11.8 

8.6 

46.1 

41.9 

-0.1 

Associated  with 
Thunderstorms 

11.4 

11*7 

40.4 

42.6 

+0.2 

*In  minutes. 


It  is  evident  that,  although  their  origins  are  different,  there  is  very 
considerable  similarity  in  the  characteristics  of  the  two  types  of  SEA.  Slight 
distinctions  can  be  identified  which  mainly  represent  the  correlation  between  the 
phasing  of  the  phenomena  on  the  two  different  frequencies;  this  connection  is 
much  closer  for  the  SEAs  due  to  thunderstorms  than  it  is  for  those  associated 
with  SIDs.  For  instance,  although  the  means  in  the  final  column  of  Table  2.2  are 
nearly  identical  and  represent  in  each  instance  an  almost  symmetrical  distri¬ 
bution  of  positive  and  negative  values  about  zero,  the  average  of  the  individual 
differences  (irrespective  of  sign)  between  the  onset  times  at  10  and  27  kc  is 
significantly  divergent  as  between  the  two  types  of  SE.^.  This  and  other  similar 
points  are  illustrated  in  Table  2.3. 
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TABLK  2. a 


COMPARISON  OF  PHASING  OF  BEHAVIOR  AT  10  AND  27  KC 
FOR  THE  TWO  TYPES  OF  SEA 


Type  of  SEA 

Mean  Difference* in 

Onset  Times  at 

10  and  27  kc** 

Mean  Difference*  in 

Rise  Times  at 

10  and  27  kc** 

Mean  Difference*  | 
in  Dui’ation  at 

10  and  27  kc** 

Associated 

with  SID 

4.4 

4.8 

17.4 

Associated  with 

Thunderstorms 

1.5 

2.2 

4.6 

*111  minutes. 

** Ignoring  sign. 


Table  2.3  indicates  that  if  there  are  appreciable  differences  at  10 
and  27  kc  in  tiie  time  phasing  of  onset,  maximum,  and  end,  then  the  SEA  concerned 
is  probably  associated  witli  an  SID,  However,  the  most  reliable  indicator  that  tlie 
SEA  is  of  solar  origin  remains  the  reciprocal  behavior  of  an  increase  in  noise  at 
27  kc  and  a  decrease  at  10  kc. 

2.3.5  Various  Aspects  of  SEA  Behavior 

The  foi'ms  of  SEAs  as  a  function  of  time  have  been  subdivided  by 
Kamada  (1960).  Kamada's  subdivision  is  somewhat  complicated  but  can  be 
simplified  into  two  classes,  A  and  B.  In  Class  A  the  rise  time  of  the  SEA  is 
much  smaller  than  tlic  decay  lime  and  the  maximum  is  a  sharp  peak.  For  Class 
B  the  rise  is  comparable  witli  that  for  Class  A,  but  is  followed  by  an  appreciable 
period  where  the  changes  from  the  maximum  value  arc  not  very  pronounced, 
until  ultimately  there  is  a  sudden  decay  to  Uio  origin;  this  takes  much  the  same 
time  as  the  rise.  Class  A  changes  usually  have  a  sawtooth  shape:  ty'pically  the 
Class  B  variations  appear  as  invci'ted  [is  or  Ws. 

Will'll  the  resniis  al  10,  2  1.  and  27  kc  are  divided  into  Classes  A  and 
B,  there  arc  sonic  iiileresi iiig.  re.’diires.  r;ible  2.']  has  been  constructed  using  63 
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SEA  inciclL'iils  whicli,  by  liiu  Ci’iterioii  of  an  increase  in  noise  on  21  anJ  27  ke  and 
a  decrease  al  10,  are  definitely  to  be  associated  with  an  SID. 


TABLE  2.4 

CLASSIFICATION  OF  SEA  INCIDENTS  ACCORDING  TO  SHAPE 
Sawtootli  (A)  Inverted  U  (B) 


Frequency  (kc)  10 

21 

27 

Number  of  Incidents 

-B 

B 

B 

13 

-B 

B 

A 

6 

-B 

A 

B 

0 

-B 

A 

A 

21 

-A 

B 

B 

6 

-A 

B 

A 

0 

-A 

A 

B 

0 

-A 

A 

A 

17 

Incidents  where  all  three  frequencies  have  the  same  class  are 
common;  this  is  probably  a  reflection  of  the  temporal  variation  in  tlie  incident 
flux.  The  class  on  21  kc  is  always  similar  either  to  that  on  27  kc  or  to  that  on 
10  kc;  this  shows  a  systematic  behavior  with  frequency.  Wlien  the  class  on  10  kc 
differs  from  that  on  27  kc,  usually  10  is  of  Class  B  and  27  of  Class  A;  this  is 
consistent  witli  the  idea  that  10  kc  is  affected  by  ionization  lower  in  the  iono¬ 
sphere  than  is  27  kc,  and  tlie  recovery  is  therefore,  perhaps,  more  sudden. 

Rough  rules  can  be  formulated  relating  the  magnitude  of  an  SEA  to 
the  parameters  of  the  causative  flare.  These  rules  should  be  taken  only  as  in¬ 
dicative;  obviously,  since  an  SID  changes  the  propagation  characteristic,  the  size 
of  an  SEA  will  depend  upon  the  distance  of  tlie  sources  of  atmospheric  noise. 
However,  the  following  average  relations  apply  to  the  Japanese  observations  for 
1958  [IGY  Data  (19.59)]. 
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At  27  Uc,  eiiliaiieeineiil  -  2.5  +  0.5  (S-i)  db 
at  10  Uc,  eiihaiicuiueni  =  -  |o.2  +  0.4  (S-1)  |  db 
where  S  is  the  importance  of  tlie  flare.  In  tlie  Japanese  equipment  the  bandwidth 
is  500  cps;  tlie  0-db  level  is  1  pr/m;  and  the  enliancement  is  defined  as  20  log 
/E;^,  where  E.^  is  tlie  maximum  intensity  of  the  SEA  and  E^  the  pre-SEA 
value. 

The  Japanese  1958  observations  were  also  examined  in  an  effort  to 
determine  any  relationship  between  solar  zenith  angle  and  the  size  of  an  SEA. 
Sueh  a  relationship  cannot  be  expected  to  be  pronounced.  Although  it  is  probable 
that  a  solar  flare  produces  its  maximum  ionospheric  effect  at  the  subsolar  point, 
the  measurement  of  atmospheric  noise  involves  integration  over  many  propaga¬ 
tion  patlis,  and  tlierefore  a  wide  area  and  a  considerable  range  of  zenith  angles. 
Table  2.5  shows  the  relationship  between  the  average  size  of  an  SEA  observed 
at  27  ke  and  Japanese  local  time. 


TABLE  2.5 

RELATIONSHIP  BETWEEN  SIZE  OF  SEA  OBSERVED 
AT  27  KC  AND  LOCAL  TIME 


Local  Time 

Number  of  SEA  Incidents 

Average  Magnitude  (db) 

800 

12 

2.8 

1000 

13 

3.5 

1200 

12 

3.0 

1400 

16 

2.8 

1600 

11 

2.4 

1800 

14 

2.4 

There  is  a  slight  tendency  for  the  larger  SEAs  to  occur  near  local 
noon,  that  is,  when  the  solar  zenith  angle  is  least,  but  the  trend  is  not  marked. 
Kamada's  (1960)  results  show  an  interesting  distinction  in  behavior  between  10, 
21,  and  27  ke.  This  is  illustrated  in  Table  2.6. 
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TABLE  2.G 


OCCURRENCE  OF  SEA  INCIDENTS  AS  A  FUNCTION  OF  LOCAL  TIME 


Number  of  SEA  Incidents  Observed 

21  kc 

27  kc 

0500-0700 

1 

20 

23 

0700-0900 

7 

113 

151 

0900-1100 

23 

118 

167 

1100-1300 

43 

112 

165 

1300-1500 

51 

137 

185 

1500-1700 

31 

137 

179 

1700-1900 

12 

81 

87 

1900-2100 

0 

7 

9 

Obviously  there  is  a  pronounced  maximum  for  the  10-kc  observa¬ 
tions  near  local  noon.  At  21  and  27  kc,  on  tlie  other  hand,  the  rate  of  occurrence 
of  SEAs  is  almost  uniform  throughout  the  day.  This  suggests  tliat  the  radiation 
responsible  for  the  10-kc  enliancement,  which  presumably  penetrates  to  lower 
heights  than  tliat  affecting  27  kc,  is  fairly  sensitive  to  tlie  value  of  the  solar 
zenitli  angle. 

One  further  aspect  of  the  Japanese  SEA  observations  is  worth  not¬ 
ing.  This  is  the  close  correlation  [Kamada  (I960)]  between  outbursts  of  solar 
radio  noise  at  ccntimetric  wavelength  (>  3000  Me)  and  SEA  incidents.  The 
start  of  the  solar  radio  outburst  usually  precedes  liiat  of  tlie  SEA  by  a  few  min¬ 
utes. 
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3.  BALANCE  EQUATIONS  FOR  THE  LOWER  IONOSPHERE 


3.1  THE  ION  BALANCE  EQUATIONS 

The  equations  governing'  the  production  and  disappearance  of  electrons  and 
ions  in  the  ionosphere  may  be  written  in  the  following  convenient  form  essen¬ 
tially  due  to  Crain  (1961): 


rf.V 

-  ()_  ,1/V_  /iV.V^  +  ic  i  /)}  A'-  (a) 

dt 

ds- 

-  1 V  -  (C  .  /))  Y-  -  (b)  (3.1) 

dt 

</W  d 

-  =  U  -  «V,\*  -  =  —  (A’  .  A'-)  .  (c) 

dt  dt 

In  these  equations  N,  N“,  and  N*  represent  the  densities,  respectively,  of  elec¬ 
trons,  negative  ions,  and  positive  ions,  while  Q  is  the  rate  of  production  for 
electrons  (and  for  positive  ions).  The  rate  coefficients.  A,  B,  C,  D,  and  E, 
govern  the  loss  and  gain  of  the  various  electrified  constituents.  A  represents 
the  loss  of  electrons  (gain  of  negative  ions)  by  the  attachment  of  free  electrons 
to  neutral  particles;  C  and  D  involve  the  reverse  processes,  being  detachment 
coefficients  for  the  removal  of  an  electron  from  a  negative  ion.  C  is  the  photo- 
detachment  coefficient,  while  D  represents  all  other  detachment  processes. 

B  is  the  electron-ion  recoml>ination  coefficient,  the  process  entailing  a  loss 
both  of  electrons  and  of  positive  ions.  E  is  the  ion-ion  recombination  coefficient, 
the  loss  now  Ijcing  of  Ixdh  positive  and  negative  ions. 
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It  is  usual,  in  attempting  solutions  of  the  balance  equations,  to  consider 
that  over-all  electrical  neutrality  is  always  maintained  and  that  Eq.  (3.1c) 
implies 


A'"  =  /V  +  N-.  (3.2) 

Q  and  the  rate  coefficients  are  all  functions  of  height.  In  addition,  Q  is  a 
function  of  time,  since  much  of  the  ionizing  radiation  is  of  solar  origin;  during 
an  SID  the  variation  of  Q  with  time  will  be  very  rapid.  Some  of  the  rate  coef¬ 
ficients,  for  example,  C,  arc  also  time  dependent  since  a  solar  flux  is  involved. 

3.2  THE  RATE  COEFFICIENTS 

Before  discussing  the  individual  coefficients  in  detail,  it  is  convenient  to 
tabulate  some  of  the  relevant  properties  of  the  atmosphere  for  the  height  range 
involved.  This  is  done  in  Table  3.1,  for  daytime  conditions. 

The  information  in  the  first  five  columns  is  essentially  that  given  in  the 
U.S.  Air  Force  ARDC  Handbook  of  Geophysics  (1960).  The  next  two  columns 
represent  a  compromise  between  the  ARDC  Atmosphere  and  the  survey  of 
Barth  (1961);  these  models  are  similar  in  that  both  give  a  maximum  of  atomic 
oxygen  concentrations  near  100  km.  There  is,  however,  a  significant  difference 
in  that  Barth  indicates  that  the  number  density  for  O  never  exceeds  that  of  02 
over  the  height  range  of  40  to  160  km;  the  ARDC  Handbook,  on  the  other  hand, 
shows  more  O  than  O2  above  90  km. 

Column  8,  the  NO  concentration,  is  taken  from  the  paper  of  Barth.  The 
figures  given  lie  between  those  of  the  ARDC  Handbook  and  those  quoted  by 
Nicolet  and  Aikin  (1960);  they  arc  approximately  two  orders  of  magnitude  less 
than  the  former,  and  one  order  of  magnitude  greater  than  the  latter. 

It  will  be  assumed  that  Table  3.1  aijplies  to  both  night  mid  day.  Barth  (1961) 
has  shown  that  the  main  alteration  with  the  change  to  night  is  a  rapid  fall  in 
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PllOPBllTlES  OP  'I'HP  ATMOSPHERE 


lit 

(km) 

Temp 

(“K) 

Pressure 

Number  Density  cm“ 

3 

(mierobar) 

Total 

N, 

O2 

0 

NO 

40 

261 

3000 

8.32  X  10^^ 

6.50  X  lO’-^ 

1.74  X  lOlf^ 

2.5  X  10^° 

2.5  X  lo"^ 

45 

276 

1600 

4.20  X  10^^ 

3.28  X  10l<5 

8.81  X  lO’-^ 

4.0  X  lOlO 

2.0  X  lO"^ 

50 

283 

879 

2.25  X  10^^ 

1.76  X  lo’^*^ 

4.72  X lO^^ 

5.0  X  10^® 

1.6  X  lO"^ 

55 

276 

484 

1.27  X  10^® 

9.91  X  lO^^ 

2.66  X  lO^S 

7.5  X  lO^O 

1.3  X  lo"^ 

60 

254 

257 

7.33  X  lOl^ 

5.72  X  10^^ 

1.54  X  10I5 

8.2  X  lOlO 

7.9  X  10® 

65 

232 

129 

4.02  X  IOI5 

3.14  X  IOI5 

8.42  X  10l4 

1.2  X  lO^l 

5.0  X  10^ 

70 

210 

60.2 

2.08  X  10I5 

1.70  X  10^^ 

4.37  X  I0I4 

1.6  X  10^1 

3.2  X  10® 

75 

188 

25.9 

1.00  X  10^^ 

7.81  X  lo’-'^ 

2.10  X  lOl"^ 

2.5  X  10^^ 

1.6  X  106 

80 

166 

10.1 

4.41  X  lO^"^ 

3.45  X  10^“^ 

9.26  X 10^2 

3.5  X  loll 

6.3  X  10^ 

85 

166 

3.69 

1.61  X  lOl'i 

1.26  X  lo’-'^ 

3.37  X  10l3 

8.1  X  lOll 

2.0  X  10^ 

90 

166 

1.35 

5.92  X  10l3 

4.64  X  10l3 

1.14  X  I0I2 

1.0  X  IOI2 

7.9  X  10^ 

95 

180 

0.512 

2.06  X  10l3 

1.61  X  10l3 

3.13  X  IOI2 

1.2  X  IOI2 

6.3  X  lO"! 

100 

199 

0.214 

7.78  X  10^2 

6.08  X  10^2 

1.00  X  10^2 

6.5  X  loll 

7.9  X  lO'l 

the  concentration  of  atomic  oxygen  at  heights  of  70  km  and  below.  This  is  not 
of  great  significance  in  the  calculation  of  rate  coefficients,  but  has  been 
included  as  a  slight  modification  in  Sec.  3.2.5. 

3.2.1  The  Attachment  Coefficient  A 

Negative  ions  are  readily  formed  in  the  atmosphere  by  the  attach¬ 
ment  of  electrons  to  atomic  (O)  or  molecular  oxygen  (O2  ).  It  does  not  appear, 
however,  that  any  of  the  other  common  atmospheric  gases  are  active  in  the 
formation  of  negative  ions.  Three  attachment  processes  occur: 
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Three -body 

e  1  0,  1 

\  . 

o,"  1  \*  .  .  . 

.  .  .  (a) 

Radiative 

e  +  Oj 

{In' . 

.  .  .  (b) 

Dissociative 

<•  +  O2 

.  0" 

.  0  . 

...  (c) 

In  Eq.  (3.3)  e  represents  an  electron;  hi/  is  a  quantum  of  radiation; 

X  denotes  a  third  body,  and  X*  is  an  excited  state  of  that  body.  The  first  two 
processes  occur  both  for  atomic  and  molecular  oxygen,  but  the  dissociative 
mechanism  is,  of  course,  confined  to  the  molecular  case. 

When  electrons  are  formed,  the  spread  in  energies  is  wide.  How¬ 
ever,  there  is  a  conversion  into  thermal  energies  of  a  few  hundredths  of  an 
electron-volt  within  a  time  of  less  than  10"®  sec.  Attachment  occurs  sub¬ 
sequently;  this  implies  that  the  process  of  dissociative  attachment,  which,  at 
the  altitudes  concerned,  is  only  comparable  with  the  three-body  process  for 
electron  energies  exceeding  some  0.5  ev,  can  be  neglected. 

The  coefficient  of  radiative  attachment  for  atomic  oxygen  has  been 
given  (  Branscomb  et  al. (1958)1  as  (1.3)  x  10"^^  cm®/sec  at  temperatures  between 
200°  and  400°  K.  The  corresponding  coefficient  for  molecular  o.xygcn  is  very 
luicertain;  it  is  believed  to  be  between  10“^®  and  10"'®. 

Biondi  (1961)  has  discussed  the  three-body  attachment  coefficient 
for  molecular  oxygen.  His  results  show  that  molecular  oxygen  is  about  thirty 
times  as  effective  as  N2  in  stabilizing  the  three-body  process.  The  dependence 
on  electron  energy  is  such  that  the  attachment  reaches  a  maximum  for  energies 
of  about  0.10  ev.  This  is  considerably  higher  than  thermal  energies  in  the 
lower  ionosphere,  so  that  electrons  initially  formed  with  energies  exceeding 
0.10  ev  have  to  pass  the  attachment  maximum  before  becoming  thermalized. 
Many  will  not  succeed  in  doing  this,  and  consequently  the  effective  attaclimcnt 
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coefficient  may  Ijc  expected  to  be  somewhat  higher  than  the  value  corresponding 
to  thermalization.  Biondi's  results  suggest  that  3  x  10'^°  is  a  not  unreasonable 
maximum  value  for  an  effective  attachment  coefficient.  Biondi  has  also  investi¬ 
gated  the  dependence  of  the  attachment  coefficient  upon  gas  temperature;  the 
results  are  somewhat  uncertain  but  indicate  that  in  the  D  region  the  coefficient 
becomes  less  as  the  temperature  decreases. 

Three-body  processes  involving  atomic  oxygen  occur  but  are  not 
significant  in  the  lower  ionosphere.  Where  atomic  oxygen  is  present  in  ap¬ 
preciable  proportions,  that  is,  above  about  90  km,  radiative  attachment  is  the 
dominant  process.  The  three-body  reaction  is  ineffective  both  because  of  the 
pressure  dependence  and  because  the  relative  concentration  of  O2,  the  most 
efficient  third  body,  is  diminished;  an  atom  such  as  O  is  more  limited  in  its 
degrees  of  freedom  than  a  molecule,  and  consequently  cannot  accept  the  excess 
energy  involved  in  Eq.  (3.3a)  as  easily. 

Summarizing,  the  attachment  coefficient.  A,  may  be  written  as 


.1  I  .3  fnl  X  10“'^  +  [0,  I  X  10“'’ 


/TjJOoX 

/ 


(oj 


[^l\3x  in-^o 
~30  j 


1 


(3.4) 


In  Eq.  (3.4)  the  symbol  [  1  represents  the  number  density  per  cm^ 
of  the  constituent  enclosed  by  the  symbol;  T  is  the  temperature  in  °K  over  the 
range  (I BO-SOO^K). 

3.2.2  The  Electron-Ion  Recombination  Coefficient  B 

Three  processes  may  be  identified:  radiative,  three-Iiody,  and 
dissociative.  In  reaction  form  thej'  may  lie  written  as 


:?.3 


Radiative 

1/'  .  <•  .  1/  ,  hi.' 

(a)  1 

Three -body 

\r  >  c  1  X  .  1/  ,  V* 

(1^)  ^ 

Dissociative 

!/*  w  .  1*  .  .D*- 

(e)  \ 

for  the  case  of  a  molecular  positive  ion  M\  A  denotes  an  atom.  The  coefficient 
for  Reaction  (3.5a)  is  about  10“'^  .  Massey  and  Burhop  (195G)  give  the  coefficient 
for  (3.5b)  as  being  2  x  10"’  at  ground  level:  since  the  three-body  process  is 
directly  proportional  to  pressure,  it  follows  that  the  three-body  coefficient 
must  be  less  than  10"®  at  heights  exceeding  40  km.  Most  estimates  for  the  dis¬ 
sociative  coefficient  arc  of  the  order  of  10"’ to  10'®;  hence  this  is  the  significant 
process  in  the  lower  ionosphere.  The  suggestion  of  Mitra  (1952)  that  the  three- 
body  process  is  dominant  below  65  km  docs  not  seem  to  be  correct. 

The  dissociative  recombination  of  NJ  has  been  examined  by  Bialecke 
and  Dougal  (1958)  and  by  Fairc  and  Champion  (1959).  Their  results  arc  in  good 
agreement  as  to  order  of  magnitude,  but  differ  in  the  form  of  dependence  on 
pressure  and  temperature.  Faire  and  Champion  find  a  value  of  4  x  10"’  at  a 
temperature  of  400°K,  and  indicate  that  the  coefficient  varies  inversely  as  the 
square  root  of  temperature,  but  is  independent  of  pressure  at  heiglits  above 
40  km;  these  results  arc  consistent  with  the  theoretical  w'ork  of  Bates  (1950). 

On  the  other  hand,  Bialecke  .and  Dougal  suggest  a  marked  and  complicated 
dependency  on  both  pressure  and  temperature. 

Biondi  and  Brown  (1949)  and  Holt  (1959)  have  examined  the  dis¬ 
sociative  recombination  process  for  Ov  .  In  eaeli  jiaper  a  value  of  about 
3  x  10"’  is  quoted. 

Ilow'ever,  it  seems  that  r.ate  coefficients  for  O^  and  nJ  have  lit'.le 
direct  application  to  the  lower  ionosphere.  Altliough  Ov  and  Nj  are  fornu'd  in 
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signiricant  quantities,  conversion  to  NO^  occurs  very  rapidly  by  I’oactions  such 


as 


■Vj"  4  0,^.  (>2^  I  ,V,  .  ,V0^  +  VO 
V,’  ,  0  .  VO^  I- 
0,/  H  ,V  .  VO^  +  0. 

Direct  measurements  [Johnson  et  al.  (1958)]  show  that  NO*  is  certainly  the 
dominant  ion  below  100  km  at  night,  and  probably  also  by  day.  The  experiments 
of  Istomin  (1961  a,b)  indicate  that  NO*  is  always  the  most  plentiful  positive  ion. 

Unfortunately  the  rate  coefficient  for  the  dissociative  recombination 
of  NO*  has  not  been  measured.  Nicolet  and  Aikin  (1960)  have  suggested  that  the 
coefficient  is  an  order  of  magnitude  less  than  the  value  for  Oj  basing  this  con¬ 
clusion  upon  the  relatively  long  lifetime  of  the  NO*  ion.  This  seems  a  reason¬ 
able  argument.  Accordingly,  B  will  be  taken  as  given  by 


/i  =  3  X  10“" 


1  '2 


(3.6) 


the  temperature  and  pressure  dependence  given  by  Faire  and  Champion  (1959) 
for  N2  being  adopted. 

Equation  (3.6)  may,  as  is  evident  above,  be  appreciably  in  error. 
An  upper  limit  for  B  is  probably  3  x  10”^,  corresponding  to  the  dissociative 
recombination  of  C^,  and  it  is  interesting  to  note  that  this  value  was  deduced  by 
Whitten  and  Poppoff  (1961)  from  ionospheric  observations.  A  lower  limit  can¬ 
not  be  assigned  at  pi’esent. 


3.2.3  The  C  Term  Photocletachment 


The  photodetachment  reactions  are  essentially 


0”  +  hi'  •  0  I  f 


Oj"  ♦-  111/  -»  Oj  +  e. 


Burch,  Smith,  and  Branscomb  (1958)  have  determined  the  coefficient  for  O2  as 
0.44;  this  process  involves  an  integration  over  the  whole  solar  spectrum. 
Although  a  value  of  0.15  ev  was  taken  in  this  work  for  the  electron  affinity, 

Ea  ,  of  the  O2  ion,  and  this  is  now  believed  to  be  low,  it  is  not  thought  that  a 
revision  of  the  calculation  with  =  0.40  --  a  more  likely  value  --  would  make 
any  significant  difference.  This  is  because  1.0  ev  corresponds  to  a  quantum  of 
wavelength  well  in  the  infrared;  since  the  peak  of  solar  radiant  energy  is  in  the 
visible,  variation  of  the  lower  limit  of  the  photodetachment  integral  over  a  wide 
range  in  the  infrared  should  have  relatively  little  effect  on  the  value  deduced 
for  C. 

Photodetachment  of  O*  can  also  occur,  the  coefficient  according  to 
Branscomb  et  al.  (1958)  being  1.4.  Thus  the  photodetachment  coefficient  is 
dependent  upon  the  proportions  of  the  negative  ions  in  atomic  and  in  molecular 
form.  Also,  since  solar  radiation  is  involved,  there  is  a  dependence  on  solar 
zenith  angle,  Y,  which  may  be  taken  after  Mitra  (1952)  as  proportional  to  cosY. 
Combining  these  effects,  C  may  be  written  as 


c 


(0-4tl0-!4  M[0-1\  _ 

- - - (tn  srr  i 

\  10 -]  4  to-]  / 


(3.7) 


and  the  proportions  of  |  O2]  and  ( 0']  may  be  considered  as  defined  by  their 
relative  rates  of  formation  by  the  attachment  processes. 
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3.2.4  The  Collisional  Detaehinciit  Coefficient  D^. 


At  heights  below  about  90  km,  the  (Tj  ion  is  the  prevalent  negative 
ion;  between  90  and  100  km,  O”  ions  arc  important,  and  indeed  in  this  region 
the  attachment  coefficient.  A,  is  dominated  by  the  radiative  attachment  process 
for  atomic  oxygen.  Collisional  detachment  is  represented  by  the  reactions 


+  1/  -•  Oj  +  ,1/  +  e 


0-  +  1/  0  +  1/  +  c. 


The  molecule  M  must  possess  sufficient  kinetic  energy  to  overcome  the  detach¬ 
ment  energy  or  electron  affinity,  E^,  of  the  ion  concerned. 

All  estimates  of  the  appropriate  electron  affinities  exceed  0.15  cv: 
that  is,  the  energy  required  to  detach  the  electron  is  considerably  larger  than 
the  mean  molecular  energy  at  temperatures  of  150°-300°K.  Thus  the  fraction 
of  the  total  number  of  molecules  that  have  the  necessary  detaching  energy  can 
be  derived  from  kinetic  theory,  using  the  approximation  for  speeds  much  greater 
than  the  most  probable  molecular  velocity.  Tlie  fraction  is  given  by 
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where  k  is  Bolzmann's  constant  and  T  the  temperature;  the  exponential  term 
is  Ijy  far  the  most  significant. 

Regarding  the  detachment  of  O2  ,  Bailey  (1959)  considered  that  all 
molecules  in  the  atmosphere  were  capable  of  producing  detachment,  but  Phelps 
and  Pack  (19G1)  find  that  molecular  oxygen  is  the  only  effective  detaching  agent. 
When  this  is  taken  into  account  and  it  is  recognized  that  the  detachment 


coefficient  is  proportional  to  Uie  number  density  and  that  the  ratio  of  Ihe  ionic 
to  tlie  molecular  free  path  is  proportional  to  the  sciuare  root  of  temperaliirc, 
the  following  exj^ression  for  Dc  for  (V.  ions  is  obtained; 


1  '2 


in® 


(3.8) 


where  P  is  the  pressure  in  microbars. 

The  value  of  remains  a  subject  of  controversy.  Twenty  years  ago  1.0  cv 
was  an  accepted  value  [Bates  and  Massey  (1943)].  More  recently,  Smith,  Burch, 
and  Branscomb  (1958)  suggested  an  affinity  of  as  low  as  0.15  ev,  but  subsequently 
there  has  been  a  steady  reaction  in  favor  of  higher  values  as  exemplified  by  the 
0.25  ev  of  Bailey  and  Branscomb  (1960),  the  0.46  ev  of  Phelps  and  Pack  (1961), 
and  the  0.74  ev  of  Jortner  and  Sokolov  (1961).  The  theoretical  work  of  Mullikcn 
(1961)  perhaps  indicates  the  start  of  a  second  reaction  towards  lower  values 
again;  Mulliken  estimates  that  lies  between  0.15  and  0.5  ev.  As  Crain  (1961) 
has  emphasized,  because  of  the  exponential  factor  a  slight  change  in  E^  implies 
a  very  large  alteration  in  Dc .  (The  same  is  true  of  changes  in  T  and  this  has 
an  interesting  consequence,  upon  the  structure  of  the  lower  ionosphere.)  Ac¬ 
cordingly,  in  Table  3.2  values  of  Dc  as  calculated  from  Eq.  (3.8)  are  given  for 
three  electron  affinities,  0.15,  0.40,  and  0.65  ev. 

As  has  been  indicated,  most  of  the  negative  ions  above  90  km  arc 
of  atomic  oxygen.  The  electron  affinity  for  the  O”  ion  is  about  1.5  ev. 

[Branscomb  et  al.  (1958)].  Thus  collisional  detachment  of  0~  is  a  most  inef¬ 
ficient  process  and  may  be  entirely  neglected.  Since  only  O2  ions  can  be 
detached,  a  correction  of  Dc  is  necessary  for  the  larger  heights  in  Tabic  3.2. 

This  is  done  approximately  by  multiplying  the  deduced  Dc  by  a  factor  of  pro¬ 
portionality  given  by  A(02]/A(02  ]  +  AjO”]  where  A|  ]  is  the  attachment  rate, 
that  is,  the  rate  of  formation,  of  the  ionic  constituent  indicatc'd  in  the  bracket. 


I'AllLH  3.2 

DFTACI! M FNT  COEFFICIE NTS 


Corrected 

1 

D, 

h(km) 

0.15 

0.40 

0.05 

(Ea  --  0.40) 

40 

2.3.10^ 

3.2 

4.9.10"^ 

3.2 

45 

1.9.10'^ 

4.4 

-4 

1.9.10 

4.4 

50 

1.0.10^ 

2.3 

-4 

1.0.10 

2.3 

55 

4 

5.7.10 

1.3 

5.7.10''^ 

1.3 

00 

2.1.10^^ 

2.9.10’^ 

4.4.10'® 

2.9.10'’’ 

05 

0.1.10^ 

2.2.10'^ 

1.3.10'^ 

2.2.10'^ 

Day  Night 

70 

1.5.10^ 

1.0.10^ 

o 

1.3.10' 

1.6.10'® 

1.6.10'®  1.6.10'^ 

75 

2.5.10^ 

4.0.10'^ 

4.5,10'^^ 

4.0.10'® 

2.5.10'®  8.1.10''’ 

80 

2.4.10 

3.3.10*^ 

5.3.10'^® 

-7 

3.1.10 

3.5.10'® 

85 

8.9 

.7 

1.2.10 

1.9.10'^® 

7.8.10'® 

8.1.10'® 

90 

3.3 

4.1.10'® 

6.9.10'^® 

q 

7.4.10' 

1.0.10'^ 

95 

2.1 

8.0.10*® 

9.3.10 

q 

1.2.10“ 

1.2.10*^ 

100 

2.1 

5.1.10'^ 

4.2.10'^® 

1.6.10'® 

6.5.10'® 

This  correction  is  applied  for  an  electron  affinity  of  0.40  cv  in  the  fifth  column 
of  Table  3.2. 

3.2..')  Associative  Detachment  (D, ) 

This  process  is  represented  by  reactions  such  as 
oT  .  0  .  n^  .  r 

and 


n"  .  0  .  0^  4  4-. 
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Estimates  of  the  rate  constants  eoncernecl  range  over  five  orders 
of  magnitude.  Whitten  and  Poppoff  (19G2)  have  shown  that  the  selection  of  a 
value  of  about  10"^"*  cm^  see“*  would  enable  some  aspects  of  the  ionospheric 
observations  of  Bailey  (1959),  during  polar  eap  absorption  events,  to  be 
interpreted  without  the  neeessity  of  postulating  significant  collisional  detach¬ 
ment  and  an  electron  affinity  of  as  low  as  0.25  ev  for  the  Og'ion.  The  laboratory 
value  of  about  0.40  ev  would  then  be  applicable. 

Although  the  uneertainties  are  very  considerable,  it  seems  best, 
at  tlie  present  juncture,  to  adopt  the  rate  constant  of  suggested  by  Whitten 
and  Poppoff  (1962).  Multiplieation  by  tlie  coneentration  of  atomic  oxygen  gives 
the  coefficients  listed  in  the  last  column  of  Tabic  3.2;  a  distinction  between  day 
and  night  is  included  for  heights  of  70  and  75  km.  At  80  km  and  above,  D  ,  >>  D, 
for  Ea  =  0.40  ev;  below  70  km  D,  exceeds  D,,  . 

3.2.G  The  E  Term  Ionic  Recombination 

There  are  two  significant  processes: 

two-lx)dy,  e.g.,  0^'  ^  o,"  -  o/  +  o^“ 

and 

lhree-lK)dy,  e.g.,  o^‘  .  .  u  >0/  <  ^ 

Values  of  between  10"“  and  3  x  10”®  have  been  quoted  for  the  two-body  process. 
Opinion  in  Britain  [Atjererombie  (1962)1  is  in  favor  of  a  figure  even  higher  than 
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3  X  10”®  ruid  Hull  in  the  U.S.A.  |Knapp  (19(i2)j  is  tending  in  the  same  direction. 
Accordingly,  the  value  of  3  x  10”®  will  be  adopted. 

The  three -body  coefficient  has  a  pronounced  dependency  upon  pres¬ 
sure  and  temperature.  It  decreases  as  the  pressure  becomes  less,  probably 
leveling  off  at  a  value  of  3  x  10”**  [Loeb  (1955)].  With  the  pressure,  P,  in 
microbars  tlie  three -body  eocfficient  may  be  written  as 

5  X  10“' 2  .  /’  .  '  3  X  10-". 


In  this  expression  the  temperature  dejxjndency  is  different,  by  a  factor  of 
from  that  given  by  Ratcliffe  and  Weekes  (1958),  The  difference  arises  because 
in  the  original  Thomson  theory  the  ratio  of  the  ionie  to  the  molecular  free  path 
is  involved;  this  is  not  independent  of  temperature  but  varies  according  to  T^^. 

Summarizing, 

/  \ 

i:  -  r,  X  10"'* .  /’  •  +  3  X  lo""  +  a  X  lo*"  (3.9) 


3.2.7  Summary 

The  rate  coefficients  arc  listed  in  Table  3.3  and  plotted  on  Fig.  3.1. 
D  represents  the  sum  of  D,  and  D., .  Crain  (19G1)  has  given  the  best  recent 
survey  of  the  data  on  rate  coeffieients;  his  work  is  especially  notable  for  its 
emphasis  upon  the  great  uncertainties  of  the  subject,  and  is  in  marked  contrast 
to  the  tendency  of  many  researchers  to  ascribe  far  more  precision  to  selected 
values  than  is  justifiable.  It  cannot  be  overemphasized  that,  although  Fig.  3.1 
is  believed  to  represent  the  optimum  development  of  the  information  presently 
(August  19G2)  available,  it  is  probable  that  within  a  year  or  two  some  of  the 
quoted  values  will  have  been  changed  by  an  order  of  magnitude! 
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FIG.  3.1  RATE  COEFFICIENTS  vs.  HEIGHT 


It  is  of  some  interest  to  examine  the  differenees  between  Fig.  3.1 
and  the  data  given  by  Crain.  As  regards  the  attaehment  eoeffieient,  A,  the  two 
sets  of  results  are  in  exeellent  agreement  exeept  for  heights  above  about  85  km 
where  the  influenee  of  atomic  oxygen  (not  considered  by  Crain)  begins  to  become 
significant.  The  B  value  is  about  40  percent  of  that  given  by  Crain.  The  results 
for  the  C  coefficient  are  identical  below  80  km,  but  above  that  level  the  effec¬ 
tive  value  for  C  is  greater  than  that  of  Crain  because  of  the  presence  of  O" 
ions.  Between  40  and  70  km  the  coefficient  D  is  primarily  due  to  collisional 
detachment;  the  agreement  with  Crain  is  reasonable,  but  the  variation  on 
Fig.  3.1  is  greater  because  of  the  inclusion  of  the  effect  of  temperature.  Above 
70  km  associative  detachment  is  the  principal  detachment  process  represented 
on  Fig.  3.1 ;  this  mechanism  was  not  considered  by  Crain,  so  that  there  is  an 
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appreciable  divergence  between  the  two  sets  of  data.  Finally,  the  E  coefficient 
is  about  forty  times  that  given  by  Crain;  this  rcllccts  the  movement  of  opinion 
referred  to  in  Sec.  3.2.6. 

3.3  The  Ionizing  Radiations  During  Quiet  Conditions 
3.3.1  Cosmic  Radiation 

Cosmic  radiation  is  usually  separated  into  the  primary  rays  incom¬ 
ing  from  extraterrestrial  sources  and  the  secondary  radiation  generated  by 
the  interaction  of  the  primary  rays  with  the  earth's  atmosphere.  The  incoming 
primary  radiation  is  deflected  towards  the  magnetic  poles  by  the  magnetic 
field  of  the  earth.  Only  the  most  energetic  primary  particles  can  avoid  such 
deflection  and  reach  equatorial  regions;  there  is  thus  a  pronounced  latitude  effect 
with  the  dependence  being  upon  a  latitude  defined  geomagnetic  ally  rather  than 
geographically.  Secondaries  are  produced  in  large  numbers  only  at  heights 
below  35  km.  However,  Meredith,  van  Allen,  and  Gottlieb  (1955)  have  shown 
that  the  ionizing  radiation  at  considerable  attitudes  consists  of  primaries, 
secondaries  generated  at  lower  heights  and  reflected  upwards  (the  splash 
albedo),  and  such  secondaries  re-entering  the  atmosphere  (the  re-entrant 
albedo).  For  geomagnetic  latitudes  of  less  than  55"  the  three  terms  are 
comparable  in  magnitude. 

Most  of  the  primary  cosmic -ray  particles  are  protons  and  have 
energies  exceeding  1  Bev.  For  such  high-energy  particles  the  rate  of  energy 
loss  per  centimeter  of  path,  and  therefore  the  energy  available  to  produce 
ionization  of  the  ambient  atmosphere,  is  almost  independent  of  the  particle 
energy.  Thus,  in  the  latitude  effect,  since  the  albedo  secondaries  are  also 
of  high  energy,  the  proportionality  is  between  ionization  production  and  the 
number  flux  rather  than  the  energy  flux  of  the  ionizing  radiation. 

The  influence  of  the  earth's  magnetic  field  means  that  incoming 
protons  with  energies  (in  Bev)  less  than  approximately  15  cos'^o,  where  o 
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is  a  ^'ecmai^aclic  laUludf,  ciuinot  ijenuirate  lo  the  latitude  o.  Thus  the  nuinlier 
flux  at  any  particular  latitude  could  be  evaluated  by  taking'  an  expression  such 
as  that  of  Ixaplon  et  al.  (1952)  for  the  energy  distribution  of  the  ineoming  primaries, 
making  allowanees  lor  the  splash  and  re-entrant  albedos,  and  integrating  be¬ 
tween  the  limits  of  15  cos’c?  and  infinity.  Refinements  eould  be  introduced 
into  the  simple  Stbrmer  eut-off  criterion  of  15  cos  “*6  as  indicated  by  Winckler 
et  al.  (1950). 

However,  the  above  procedure  although  complicated  is  fundamentally 
empirical,  as  there  is  no  sound  theoretical  basis  for  the  energy  spectrum  of 
the  incident  cosmic  radiation.  Since  an  empirical  approach  is  inevitably 
necessary,  it  seems  easiest  to  apply  it  in  the  following  way.  Measurements 
on  the  latitudinal  dependence  of  the  vertical  flux  of  particles  near  the  top  of 
the  atmosphere  have  been  summarized  by  Curtis  (1956).  Between  geomagnetic 
latitudes  of  0°  and  60°  the  latitudinal  dependence  is  fairly  well  approximated  by 
a  simple  proportionality  with  eos"^©;  at  higher  latitudes  the  flux  is  roughly 
constant.  The  loss  of  energy  of  the  particles  by  ionization  is  taken  to  be  2.2 
X  10^  ev  per  centimeter  traveled,  at  atmospheric  pressure,  and  32  ev  is 
accepted  as  the  energy  corresponding  to  the  production  by  cosmie-ray  particles 
of  an  ion  pair  in  air.  By  using  these  figures,  the  rate  of  ionization  production, 
q  ,  in  ion  pairs  per  cm*  per  second  per  atmosphere  is  given  by 

9 

/7  2:  *^8  rn.'v  *  ^  .  (3.10) 

'  n 

Equation  (3.10)  is  in  good  agreement  with  three  values  quoted  by 
van  Allen  (1952)  as  showm  in  Table  3.4. 

TABLE  3.4 

COSMIC-RAY  IONIZATION  AS  A  FUNCTION  OF  GEOMAGNETIC  LATITUDE 


() 

3° 

38° 

51° 

(I  (\'an  Allen) 

*  O 

35 

90 

257 

iFq.  (3.10)] 

38 

9.8 

243 

The  rate  of  ionization  production,  q  ,  due  to  cosmic  rays  at  a  height  h  is 

C 

readily  obtained  from  Eq,  (3.10)  since 


(lAh,(f)  =  </ 


=  38  cos  *  0  - 


38  , 


(2.S,3)  X  10* 


where  and  are  the  number  densities  at  heights  h  and  0 . 

It  is  interesting  to  note  that  the  differential  number  speetmm,  dN, 
for  incoming  primaries  of  energy  E  is  quoted  as  dN  =  kE"“dE,  where  k  is  a 
constant  and  the  empirieal  constantahas  a  value  of  between  1.8  and  1.9.  It 
follows  that  if  the  shadow  cone  effeet  is  ignored  and  the  simple  Stdrmer 

foo 

cut-off  is  adopted,  then  =  J  kF.~°flE,  which  is  not  far 

1 5  cos ^  0 

removed  from  the  cos'^e  dependency  of  Eq.  (3.10). 

One  important  feature  of  the  cosmie-ray  latitude  dependency  is 
the  system  of  geomagnetic  eo-ordinates  to  be  selected.  Recent  information 
[Sandstrdm  (1958)]  suggests  that  the  cosmic-ray  equator  corresponds  best  with 
the  dip  equator.  Thus  the  ion  production  function  is  best  written  as 

^  {h,  /)  =  1 1)  ni^cos~^l  X  lO""*.  (3-11) 


where  I  is  the  angle  of  dip. 

Regular  temporal  changes  in  cosmic-ray  activity  are  slight,  with 
the  exception  of  the  variations  associated  with  the  solar  eyele  [Simpson  (1958)]. 

It  appears  that  many  more  low-encrgy  protons  reach  the  earth  during  solar  mini¬ 
mum  than  at  times  of  maximum.  These  protons  are  naturally  channeled  into 
high  latitudes  where  the  cyclical  effects  arc  thus  likely  to  be  most  pronounced. 
The  vari.ation  in  ion  production  near  the  poles  is  probably  ±25  percent  bctv'een 
solar  maximum  and  minimum;  for  latitudes  less  than  55°,  however,  it  is 
mucli  less. 
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ion  protliiciioii  rules  (ieriveci  Ironi  Eq.  (3.11)  are  given  in  Table  3.5. 


3 .3 .2  Lynnm  o 

The  sun  emits  intensely  at  the  wavelength  of  1215.7A,  the  first  line 
in  the  Lyman  scries  of  hydrogen.  It  was  thought  a  few  years  ago  that  the 
emission  was  extremely  variable  and  that  many  of  the  short-term  fluctuations 
in  the  ionosphere,  such  as  SIDs,  were  associated  with  rapid  changes  in  the 
intensity  of  the  Lyman  -O'  radiation.  However,  as  additional  information  by  more 
reliable  methods  has  become  available,  notably  through  the  NRL  rocket  and 
satellite  program  (Friedman  (1961)1  it  has  become  evident  that  on  a  short¬ 
term  basis  the  flux  of  Lyman  o  is  remarkably  constant,  but  the  possibility  of 
a  long-term  variation,  perhaps  by  as  much  as  an  order  of  magnitude,  in 
phase  with  the  solar  cycle,  cannot  at  present  be  eliminated.  During  the  solar 
maximum  the  flux  is  about  5  ergs  cm’^sec'^at  the  top  of  the  terrestrial 
atmosphere.  This  corresponds  to  3  x  10^^  photons  cm“^sec‘\ 

Lyman-a  radiation  is  capable  of  ionizing  NO  and  is  strongly 
absorbed  by  molecular  oxygen.  Watanabe,  Inn,  and  Zelikoff  (1953)  give  the 
absorption  cross  section  for  molecular  oxygen  as  1.04  x  10"^°cm^,  while  for 
nitric  oxide  the  ionization  cross  section  is  2  x  10"^®  cm^  and  the  absorption 
cross  section  2.4  x  10“^®cm^  (Watanabe  (1954)].  For  molecular  nitrogen,  Nj, 
the  Lyman-o;  absorption  coefficient  is  probably  less  than  6  x  10’*®cm^ 

(Ditchbum  etal.  (1954)]. 

If  the  incident  flux  of  Lyman  ais  I^  ,  then  the  flux,  Ii,,  at  any  height 
h  may  be  written  as 


(3.12) 


whore  Uj  and  n^  arc  the  absorption  cross  section  and  number  density  respec¬ 
tively,  of  the  constituent  j.  F.quation  (3.12)  includes  the  approximation  of 
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replaeinj^'  llu-  Clia))man  liinctioii  /.(.‘iiilh  V  by  sec  ,V;  this  is  reasonable 

lor  X'  :i80°.  'I'he  production,  ,  oi  ionizat'on  at  height  h  is  given  by 


‘h.  "Z  " i I,'* 

; 

t 

where  Oj  represents  the  ionization  cross  section. 

Of  the  various  tenns  hi  Eqs.  (3,12)  and  (3,13),  only  that  considering 
NOis  involved  in  the  ionization  summation.  As  regards  the  absorption  summation, 
the  absorption  is  very  slight  above  100  km,  so  that  I  „  ~  Iioo"  Below  100  km,  a 
combination  of  the  information  on  Table  3.1  and  on  the  absorption  cross  sections 
shows  that  the  absorption  is  almost  entirely  due  to  molecular  oxygen.  For  the 

1  /’ 

molecular  oxygen  constituent  I  n .  <lh  =  — -  — -  where  M_  is 

the  molecular  weight  and  the  atmospheric  pressure  at  height  h.  With  g 

taken  as  060  cm  sec  over  the  range  of  heights  concerned  and  =  5.31  x 
-22 

10  gram,  Eq.  (3.12)  can  be  rewritten  as 

'h  'oo V  •  (3-14) 


By  using  Eqs.  (3.13)  and  (3.14)  and  the  information  of  Table  3.1,  Table  3.6 
giving  ion  production  rates  due  to  Lyman  o?  may  be  calculated. 

Table  3.6  has  not  been  continued  below  about  60  km,  because 
cosmic-ray  ionization  is  dominant  at  these  lower  heights.  The  minimum  at 
95  km  follows  from  the  minimum  of  NO  indicated  by  Barth  (1961)  at  this 
height. 

3.3.3  X  Rays 

X  rays  are  capable  of  ionizing  all  constituents  of  the  atmosphere, 
and  it  is  usually  taken  that  all  the  absorbed  energy  is  converted  into  ionization. 
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The  eoniputatioii  ol'  X-iniy  ioniziitiou  is  somewhat  nioi’e  difficult  than  that  for 
Lyman  oi  because  the  radiation  is  not  monochromatic,  but  the  basic  procedure 
is  essentially  similar. 


If  1  is  the  intensity  of  the  incident  radiation,  then 

r  1 

V  f‘„,  I  where  Pm  is  the  mass  absorption  coefficient, 

A  -1 

ph  th('  density  at  height  h,  and  the  secant  approximation  to  the  Chapman  function 

r 

is  used.  Now  I  p^M  ^ - ,  so  that  the  rate,  Qh  of  ion  production  at  height  h 

\  P 

may  be  expressed  as 


-  Ph  '■*/’ 


(3.15) 


where ^  is  the  yield  in  ions  per  unit  absorption. 

The  quantities  Pm  and  /3  (to  a  lesser  extent)  are  functions  of  wave¬ 
length,  while  the  incident  solar  X-ray  flux  contains  energy  over  a  wide  range  of 
wavelengths.  In  order  to  obtain  solutions  for  qn,  three  approaches  are  possible. 
First,  elaborate  numerical  integration  can  be  performed.  Second,  attempts  can 
be  made  to  obtain  functional  representations  of  the  quantities  that  are  frequency 
dependent,  with  a  view  to  obtaining  an  expression  for  qh  in  terms  of  wavelength 
X.  Third,  approximations  can  be  included,  where  appropriate,  in  an  effort  to 
obtain  solutions  that  are  reasonably  accurate  but  still  tractable  without  the 
necessity  for  extensive  computation;  this  is  the  approach  that  will  be  adopted. 

Rajewsky  and  Lang  (1961)  have  reviewed  the  information  relating 
the  average  energy  lost  in  creating  an  ion  pair  to  the  energy  of  the  incident 
X-rays.  Tlie  trend  of  the  results  is  indicated  in  Table  3.7. 
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TABLK  3.7 


IONIZING  EFFECTIVENESS  OF  X  RAYS 


Energy  of  incident  X  rays  (kev) 

0.1 

1 

10 

100 

Corresponding  wavelength  (A) 

127 

12.7 

1.27 

0.127 

Energy  lost  in  producing  one 

ion  pair  (ev) 

50 

37.5 

34.5 

33 

The  energy  loss  per  ion  pair  varies  appreciably  with  wavelength 
for  the  softer  X  radiation  but  tends  towards  constancy  as  the  X  rays  become 
harder,  Rajewsky  and  Lang  find  from  their  experiments  that,  between  2,2 
and  0.7 A,  35  ev  is  a  good  approximation  to  the  energy  loss  per  ion  pair,  and 
this  value  will  be  taken  as  applying  over  the  whole  X-ray  spectrum.  The 
error  in  this  assumption  is  appreciable  at  wavelengths  greater  than  about  15 
A,  but  this  wavelength  range  is  mostly  absorbed  above  100  km  and  is, 
therefore,  of  minor  significance  as  regards  the  lower  ionosphere. 

A  value  indicative  of  jS  having  been  established,  it  remains  to 
consider  •  Appropriate  information  concerning  as  a  function  of  wave¬ 
length  has  been  given  by  Victoreen  (1949)  and  Henke  et  al.  (1957).  A  com¬ 
posite  representation  of  tlie  available  data  is  given  in  the  curve  of  Fig.  3.2. 
For  wavelengths  exceeding  0.6  A  and  less  than  about  20  A ,  that  is,  over  the 
range  most  important  in  the  lower  ionosphere,  the  variation  in  Fig.  3.2  is 
represented  by  the  empirical  relationship  =-  (3.1)\'^-'”  .  Below 

0.25  A ,  a  relationship  witli  similar  form  but  different  constants  applies: 
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FIG.  3.2  MASS  ABSORPTION  COEFFICIENT  FOR  X  RAYS  IN  AIR  AS  A  FUNCTION 
OF  WAVELENGTH 


If  an  incident  flux,  I^,  of  1  erg  cm"®  sec*^  is  taken,  with  a  conversion 
factor  of  35  ev,  Eq.  (3.15)  reduces  to 


,  a  ,„!0  /  y  ^h\ 

1 .8  X  10  «  p.  exp  [ - I 

\  960  / 


(3.16) 
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where  960  em  sec'“  is  taken  as  a  mean  value  of  the  aceeleration  due  to 
gravity  over  the  altitude  range  eoncemed.  Equation  (3,16)  is  plotteil  as  a 
function  of  height  for  five  values  of  V  and  seven  values  of  the  wavelength 
of  the  incoming  radiation,  in  Figs.  3.3  through  3.8.  It  is  evident  that  the  harder 
the  X  radiation  the  lower  the  height  at  which  the  maximum  production  of  ioni/.ation 
occurs.  Also,  for  a  given  wavelength,  a  change  of  zenith  angle  has  a  much  larger 
effect  at  heights  below  that  of  maximum  production  than  at  those  above. 

It  remains  to  determine  the  X-ray  flux  associated  with  a  quiet  sun. 
Poppoff  and  Whitten  (1962)  have  estimated  the  background  X-ray  flux  from  data 
obtained  by  the  Sunflare  II  Probe  of  the  Naval  Research  Laboratory.  The  latest 
available  results  from  the  NRL  scries  of  rocket  and  satellite  experiments  show 
that  near  solar  maximum  the  quiet  solar  flux  in  the  2  to  8  A  band  never  exceeds 
6  X  10“'*  erg  cm“^  sec'*  .  [Kreplin  (1962)]  .  In  order  to  make  tlie  Poppoff- 
Whitten  spectrum  agree  with  this  observation  of  Kreplin,  it  must  be  reduced  by 
a  factor  of  about  2.5;  this  has  been  done  and  the  result  is  plotted  in  Fig.  3.9. 

In  order  to  determine  accurately  the  rate  of  ion  production  represented 
by  the  incident  flux,  a  complicated  calculation  involving  summation  over  all 
wavelengths  is  necessary.  However,  a  reasonable  approximation  may  be  made 
by  splitting  the  flux  into  different  wavelength  ranges  and  then  applying  the 
information  of  Fig.  3.3  to  3.7.  A  convenient  subdivision  is  shown  on  Fig.  3.9 
and  also  given  in  Table  3.8.  By  using  this  subdivision,  the  rate  of  production  of 
ionization  for  the  quiet  solar  X-ray  flux  has  been  calculated  for  several  zenith 
angles.  The  results  are  given  in  Table  3.9. 

3.3.4  Other  Ionizing  Agents-- Watanabe  and  Hintcregger  (1962)  have 
examined  the  ionization  due  to  solar  ultraviolet  radiation  and  soft  X  rays.  Their 
results  show  that  appreciable  ionization  can  be  produced  at  heights  of  95  to  100 
km  by  radiations  between  911  and  1027  A,  and  by  X  rays  with  wavelengths  in  the 
range  of  30  to  170  A.  At  90  km  and  below,  however,  this  ionization  is  insignif¬ 
icant  as  compared  with  that  already  considered  due  to  X  rays  harder  than  30  A  . 
Tabic  3.10  indicates  the  rate  of  production  of  ions  for  these  additional  radiations. 


54 


FIG.  3.3  ION  PRODUCTION  DUE  TO  X  RAYS  AS  A  FUNCTION  OF  HEIGHT  AND 
WAVELENGTH  (Zenith  Angle  0") 
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HE  tGMT  —  kiloHieterj 

FIG.  3.4  ION  PRODUCTION  DUE  TO  X  RAYS  AS  A  FUNCTION  OF  HEIGHT  AND 
wavelength  (Zenith  Angle  20") 
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FIG.  3,5  ION  PRODUCTION  DUE  TO  X  RAYS  AS  A  FUNCTION  OF  HEIGHT  AND 
WAVELENGTH  (Zenith  Angle  40  ) 
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FIG.  3.6  ION  PRODUCTION  DUE  TO  X  RAYS  AS  A  FUNCTION  OF  HEIGHT  AND 
WAVELENGTH  (Zenith  Angle  60  ) 
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FIG.  3.8  ION  PRODUCTION  DUE  TO  2A  X  RAYS  AS  A  FUNCTION  OF  HEIGHT  FOR 
FIVE  ZENITH  ANGLES 
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TABLE  3.8 


SUBDIVISION  OF  X-RAY  FLUX  FROM  QUIET  SUN 


Wavelength  Range 

Center  Wavelength 

(A) 

(A) 

19-31 

25 

1.4.10  ^ 

O 

1 

to 

13-19 

16 

-4 

9.0.10 

-3 

5.4.10 

7-13 

10 

-4 

4.5.10 

-3 

2.7.10 

5-7 

6 

1.3.10"'^ 

-4 

2.6.10 

4 

-5 

3.10 

0.10~''’ 

2 

10-^' 

2.10~‘’ 

The  "average"  meteor  appears  at  a  height  of  about  110  km  and  dis¬ 
appears  at  80  km.  Over  this  range  of  altitude  the  mean  rate  of  ion  production 
for  all  the  sporadic  meteors  occurring  globally  is  about  10  cm'  *  sec  *' 

(Lovell  (1957)] .  Taking  into  account  the  height  distribution  of  the  ionization 
yields  Table  3.11.  It  must  be  emphasized  that  meteor  ionization  is  not 
likely  to  be  uniform  below  100  km;  it  is  improbable  that  the  diffusion  of  the 
individual  trails  will  be  sufficient  for  them  to  merge  into  an  isotropic  wliolc. 
Accordingly,  it  is  somewhat  dubious  whether  meteor  ionization  should  ever  be 
regarded  as  a  continuous  source.  The  average  flu,\  of  shower  meteors  is  only 
about  20  percent  of  that  associated  with  sporadic  meteors:  however,  during  the 
largest  showers  the  ion  production  rale  may  rise  to  as  much  as  a  tliousand 
times  the  rates  given  in  Table  3.10, 
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tablp:  3.10 


RATH  OF  rON  PRODUCTION  (cm'^scc'O  BY  ADDITIONAL  RADIATIONS 


Solar  Zenith  Angle  y 

Height 

(km) 

95 

100 

0 

2.4x10 

1.3x10'^ 

20 

2.1x10 

1.2x10'^ 

40 

1.4x10 

3 

1.0x10 

60 

6.7 

2 

4.0x10 

80 

7.0xl0”^ 

4.0x10 

3.3.5  Summary — Figure  3.10  summarizes  the  ionizing  contributions  of  the 
various  radiations.  It  can  be  seen  that  below  some  60  km  the  cosmic  ray  source 
is  dominant;  above  80  km  X  rays  are  most  important;  while  between  60  and  80  km 
the  contributions  due,  respectively,  to  Lyman  a,  X  rays,  and  cosmic  radiation 
are  of  comparable  size.  The  results  of  Fig.  3.10  confirm  the  conclusion  of 
Poppoff  and  Whitten  (1962)  that  X  rays  are  the  major  factor  in  the  formation  of 
the  D  region:  the  ojjinion  of  Nicolet  and  Aikin  (1960)  that  Lyman  o  is  the  domin¬ 
ant  influence  is  not  supported.  However,  it  should  be  emphasized  that  the 


TABLE  3.11 

HATE  f)F  ION  PRODUCTION  (cm"’  sec'O  BY  METEORS 


Height 

Ion  Production 

80 

0 

-4 

85 

2.5x10 

-4 

90 

1.7x10 

-4 

95 

1.0x10 

-4 

1  00 

0.5x10 

FIG.  3.9  X-RAY  FLUX  FROM  QUIET  SUN  NEAR  MAXIMUM  OF  SOLAR  CYCLE 


X-ray  flux  from  which  the  X-ray  contribution  to  Fig.  3.10  was  derived  was 
obtained  for  quiet  solar  conditions  during  a  period  not  far  removed  from  sun¬ 
spot  maximum.  If  the  quiet  X-ray  flux  follows  the  solar  cycle  while  that  for 
Lyman  a  remains  unchanged,  then  the  relative  importance  of  X  rays  and  Lyman 
G  could  be  very  different  at  sunspot  minimum  and  at  sunspot  maximum.  Any 
conclusions  on  this  point  are  still  somewhat  indefinite,  but  it  is  very  significant 
that  Kreplin,  Chubb,  and  Friedmann  (1962)  in  a  recent  paper  estimate  that  the 
X-ray  flux  in  the  important  band  of  2-8A  changes  by  a  factor  of  a  thousand  dur¬ 
ing  the  solar  cycle;  the  ultraviolet  emission,  on  the  other  hand,  seems  relatively 
stable. 

Above  about  95  km,  as  Fig.  3.10  shows,  the  radiations  forming  the  E 
layer,  namely  solar  ultraviolet  and  soft  X  rays,  are  dominant.  Below  some  70 
km,  cosmic  rays,  as  pointed  out  by  Moler  (1960),  are  very  important,  and  even 
at  greater  heights  seem  to  produce  somewhat  larger  effects  than  those  due  to 
meteors. 
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FIG.  3.10  RATE  OF  ION  PRODUCTION  AS  A  FUNCTION  OF  HEIGHT  UNDER  QUIET 
SOLAR  CONDITIONS 
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4.  SOIXITION  OF  THE  BALANCE  EQUATIONS 
FOR  QUIET  CONDITIONS 

4.1  SOLUTION  ASSUMING  BALANCE  IS  IMMEDIATELY  ATTAINED 

The  survey  of  the  rate  coefficients  made  in  Sec.  3.2  shows  that  the  two 
recombination  coefficients  B  and  E  are  approximately  equal  over  the  height 
range  of  40  to  100  km.  An  assumption  that  B  =  E  ~  4.10"®  leads  to  a  very 
considerable  simplification  of  the  equations,  and  to  convenient  solutions;  it  will 
therefore  be  adopted.  Admittedly,  the  assumption  involves  approximation,  but 
the  resulting  solutions  are  far  more  accurate  than  some  previous  machine 
computations  which  employed  values  of  the  rate  coefficients  now  demonstrably 
incorrect.  Generally,  in  this  field,  there  is  a  dangerous  tendency  to  accept  an 
elaborate  machine  solution  as  being  entirely  correct,  whereas  there  may  well 
be  uncertainties  of  an  order  of  magnitude  in  some  of  the  parameters  upon  which 
the  solution  is  based. 

Under  quiet  conditions  the  rate  of  change  of  the  concentrations  of  electrons 
and  ions  is  small.  Thus,  by  equating  the  derivatives  to  zero  and  putting 
B  =  E,  Eq.  (3.1)  simplifies  to 


n  -  ,1V  -  «\,v*  ^  {c  in  ,\"  =  n 
IV  -  {c  *-  in  -  it\-\’*  =.  0 
o  _  /jvv*  -  ^  n  _  i{{\*)^  .  0 


(4.1) 
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The  equilibrium  solutions  are 
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where  F,  the  complete  detachment  coefficient,  is  given  by  F  =  (C  +  D). 

Profiles  of  N  as  functions  of  height  and  time  may  be  derived  using  the 
information  of  the  tables  that  have  been  given  and  interpolating  where  necessary. 
Two  illustrative  examples  will  be  given. 

The  first  of  these  is  for  an  equinoctial  day  and  a  path  lying  around  the 
equator.  The  sun  will  be  overhead  at  local  noon,  and  it  is  assumed,  incorrectly 
but  for  simplicity  in  the  computations,  that  the  dip  angle  is  zero  over  the  entire 
equatorial  path.  Figure  4.1  is  a  contour  representation  of  the  electron  distribution; 
it  is,  of  course,  symmetrical.  Figure  4.2  shows  profiles  of  N  against  height  for 
the  vertical  paths  and  zenith  angles  indicated. 

The  second  example,  on  Figs.  4.3  through  4.5  is  for  the  winter  solstice 
with  noon  over  Capricorn  and  a  meridional  path.  It  is  taken  that  the  dip  angle  is 
equal  to  the  geographic  latitude.  Evidently  the  combined  effects  of  cosmic-ray 
and  solar  ionization  introduce  considerable  asymmetries,  especially  at  the  lower 
heights. 
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FIG.  4.1  EQUATORIAL  PATH  AT  EQUINOX  -  ELECTRON  DENSITY  CONTOURS  (cm"^) 
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lAT  ITUDE 

FIG.  4.3  MERIDIONAL  PATH  AT  WINTER  SOLSTICE  -  ELECTRON  DENSITY  CONTOURS  (cm-^) 
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FIG.  4.4  PROFILES  OF  ELECTRON  DENSITY  vs.  HEIGHT  FROM  FIG.  4.3 
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4.2  DECAY  OF  IONIZATION  AFTER  SUNSET 

The  assumption  that  equilibrium  is  reached  seems  reasonably  valid  for 
daytime  hours.  Examination  of  the  rate  coefficients  shows  that  after  any  day¬ 
time  change  equilibrium  conditions  ai’e  always  re-approached  fairly  rapidly. 
Direct  exiocrimental  evidence  of  D-layer  effects,  such  as  absorption  [Bibl  ^ 
al.  (1961)]  and  VLF  phase  changes  at  near  vertical  incidence  [  Bracewell  et 
(1951)]  ,  indicates  that  a  dependence  on  solar  zenith  angle  is  followed  with  little 
time-lag.  This  again  confirms  the  concept  that  the  quiet  daytime  D  region  is 
never  far  removed  from  an  equilibrium  state. 

However,  equilibrium  is  far  from  being  attained  at  night,  particularly  at 
the  greater  altitudes.  The  ionization  remaining  after  sunset  decays  slowly  and 
the  low  equilibrium  values  are  hardly  approached  before  sunrise  occurs.  The 
contours  of  Figs.  4.1  and  4.3  have  therefore  been  modified  ignoring  solar 
radiation  for  X  >  80°.  The  values  of  the  electron  density  have  been  obtained  by 
considering  the  appropriate  nighttime  values  of  the  parameters  A,  B,  F,  and  Q 
for  any  position,  and  the  time  since  X  =  80°  at  that  position.  The  decay  of 
ionization  is  essentially  described  by  Eq.  (3.1c)  since  the  balance  of  attachment 
and  detachment  processes  is  achieved  rapidly  after  sunset  and  the  ratio  N  /N 
then  remains  constant  and  equal  to  A/F.  This  condition  and  the  condition  of 
electrical  neutrality  N  +  N  =  are  sufficient  to  specify  N,  once  the  behavior 
of  N  is  known.  Equation  (3.1c)  can  be  solved  exactly,  giving 

'V  [2.'V;^A„(W>^.xp(-2VWA)] 

'  “  ZN*  +  N^(]-  €xp{-2 ^/JQ  t)) 

4-  + 

where  N  q  is  here  the  nighttime  equilibrium  value  of  ionization  and  (N^^  +  N  ^^ ) 
is  the  equilibrium  value  existing  at  time  t  =  o,  when  X  =  80°.  The  modified 
distributions  for  N  are  shown  on  Figs.  4.6  and  4.8  together  with  cross  sections 
on  Figs.  4.7  and  4.9. 
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ELECTRON  DENSITY  -  cm" 


FIG.  4.7  PROFILES  OF  ELECTRON  DENSITY  vs.  HEIGHT  FROM  FIG.  4.6 
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L*TITUDE 

FIG.  4.8  MERIDIONAL  PATH  AT  WINTER  SOLSTICE  CONSIDERING  NONEQUILIBRIUM 
AT  NIGHT  -  ELECTRON  DENSITY  CONTOURS  (cm’^) 
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ELECTS 


FIG.  4.9  PROFILES  OF 


DENSITY  vs.  HEIGHT  FROM  FIG.  4.8 


4.3  COMPARISON  WITH  EXPERIMENTAL  OBSERVATIONS 


Several  theoretical  models  of  electron  density  profiles  in  the  lower 
ionosphere  have  been  given,  but  actual  ejqjerimental  observations  are  very  scanty. 
Profiles  may  be  deduced  by  radio  probing  from  the  ground  and  examination  of 
the  partial  reflections  returned  from  the  D  region;  this  approach  has  been 
employed  recently  by  Titheridge  (19G2).  The  direct  method  of  measuring 
electron  densities  using  rocket  probes  is  difficult  for  heights  below  100  km, 
since  the  density  of  the  neutral  particles  is  so  much  greater  than  that  of  the 
electrons.  However  Smith  (1962)  has  reported  three  successful  rocket  flights. 
Some  of  the  experimental  results  are  given  on  Fig.  4.10  with  the  appropriate 
theoretical  profiles  deduced  from  Fig.  4.8  for  comparison.  The  agreement  is 
reasonably  good.  It  would  be  surprising  if  it  were  exact,  since  the  theoretical 


FIG.  4.10  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  CURVES 
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work  lui.s  comiJletoly  i^noreil  Iho  vertical  transport  ot  electrons  by  diffusion  or 
other  processes.  Perliaps  the  major  differeiiee  between  the  theoretical  and  the 
e.xi)eri mental  curves  is  the  tendeiiey  of  the  former  to  give  greater  electron 
densities  than  the  latter  at  heights  near  80  km;  this  is  true  both  by  night  and  by 
day. 

4.4  RELATION  OF  ELECTRON  DENSITY  PROFILES  TO  VLF  PROPAGATION 

The  propagation  characteristics  for  VLF  radio  waves  in  the  quasi-waveguide 
formed  by  the  surface  of  the  earth  and  the  lower  ionosphere  are  chiefly  deter¬ 
mined  by  the  variation  in  ionospheric  conductivity  with  height.  For  this  reason 
it  is  desirable  to  convert  the  contours  of  electron  density  into  contours  of  the 
parameter  .  This  parameter  is  directly  related  to  ionospheric  conductivity 
and  is  commonly  employed  in  the  development  of  the  waveguide  theory  of  VLF 
propagation.  The  equation 


1/2 


(4.4) 


defines  ;  in  Eq.  (4.4),  vis  the  electron  collisional  frequency,  the 
longitudinal  component  of  the  electron  gyrofrequency,  and  the  plasma 
frequency.  In  mks  units,  where  e  is  the  electronic  charge 

(1.6  X  lO”'”  c),  m  is  the  electronic  mass  (9.1  x  10"^'  kg),  8.85  x  10“'^, 
and  N'  is  the  number  of  electrons  per  m  ^  (N'  =  10  ®  N).  Hence  =  3.4  x  10®  N. 
The  quantity  of  is  given  by  (4"  x  10“^  )He/m  where  H  (ampere-turn/meter) 
is  the  strength  of  the  longitudinal  component  of  the  earth's  magnetic  field  for 
the  propagation.  The  best  information  upon  collisional  frequencies  is  in  the 
report  by  Cook  and  Lorents  (1961)  (sec  Fig.  4.11)  while  the  paper  of  Mlodnosky 
and  Hclliwcll  (1962)  is  useful  in  determining  for  meridional  propagation. 
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The  electron  densities  shown  on  I'ifiS.  4.G  and  4.8  have  been  converted  to 
contours  of  ;  the  results  are  given  on  Figs.  4.12  through  4.19.  For  the 
equatorial  case,  was  plotted  for  propagation  along  the  equator  ( -  0); 
the  meridional  case  is  given  both  for  0,  and,  more  correctly,  for  the 

instance  of  propagation  along  the  meridian  ( non-zero).  In  all  cases  the 
divergence  between  magnetic  and  geographic  north  is  neglected. 

It  is  noteworthy  that  the  contours  of  are  far  from  being  parallel  to  the 
surface  of  the  earth.  This  implies  that  in  the  waveguide  theory  of  VLF  propagation 
the  guide  must  be  considered  to  be  of  variable  width.  The  effects  are  briefly 
discussed  in  Appendix  11. 
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4.18  PROFILES  OF 


HFIGHT  FROM  FIG.  -1. 17 
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Tho  SID,  as  already  indicaled,  is  due  to  an  increase  in  X-ray  tktx,  that  is. 
an  enhancement  ol  the  flux  pjiven  by  Fij;'.  .’i.i).  The  intensity  anrl  spectral  distri¬ 
bution  of  the  X  radiation  uiifloubtodly  varies  between  different  SID  ineiilenls,  t)ul 
any  investigation  of  such  variation  must  await  the  analysis  ot  extensive  satellite 
data.  Indeed,  it  ean  hardly  be  said  that  even  the  f>enei';il  behavior  ftiiriiifr  an  SID 
has  yet  been  established. 

.'i.l  SPECTRAL  DISTRIBUTION  OF  X-RAY  FLUX 

Poppoif  and  Whitten  (1961)  have  dedueed  the  X-ray  speeti-um,  associated 
with  a  Class  2+  flare,  from  rocket  observations:  the  resulting  D-l.ayor  electron 
distribution  produces  ionospheric  absorption  that  is  in  reasonable  agreement 
with  that  experimentally  observed.  Kreplin  (1961)  states  that  rluring  a  fUire  the 
X-ray  emission  in  the  2  to  8A  band  inereases  by  a  factor  of  ten  to  one  hundi'ed: 
it  can  be  deduced  from  the  same  paper  that  over  the  wa'  clength  range  of  8  to 
20A  the  faetor  is  three  to  twenty-five,  while  in  the  band  between  41  and  60A  thi' 

SED  flux  is  probably  only  about  three  times  --  and  certainly  less  than  ten  times  -- 
as  strong  as  that  in  undisturbed  conditions.  The  satcdlite  observations  of  Kreplin. 
Chubb,  and  Friedmann  (1962)  show  that  during  a  Class  2  flare  aecomiianied  l)v  a 
short-wave  faflcoul  (also  of  Class  2)  the  X-ray  flux  in  the  2  to  8A  Ixind  is  some 
T)  X  I0“’^lo  10"^  erg  cm^  ''see. 

The  papers  of  Kreplin  and  his  colleagues  suggc'st  that  the  Whitten -Popi)off 
X-ray  spectrum  is  considerably  more  intense  th;m  that  associated  with  a  mod- 
ei'ate  flare.  In  order  to  reconcile  the  quiet-sun  and  solar-flni'c  ol)servations. 
it  is  necessary  to  reduce  the  magnitude  of  the  X-ray  Ilux  indicatc'd  l)v  Whittem 
and  Poppoff  by  a  factor  of  thirty,  whih'  retaining  the  same  stx'ctrnl  <tist  rilnit  ion. 
This  has  been  done  and  the  result,  idotted  in  i'Tg.  a. I,  will  lu'  tal-.en  as  represtmt  in 
the  X-ray  flux,  addition.al  to  lh('  (|uiet  \alue,  r('sulting  from  a  mod('rate  Hare. 


FIG.  5.1  ADDITIONAL  X-RAY  FLUX  DUE  TO  A  MODERATE  FLARE 


The  X-ray  flux  during  a  flare  rises  rapidly  to  a  maximum  and  then  decays 
more  gradually.  It  is  probable  that  the  predominant  X  radiation  steadily  becomes 
softer  as  time  increases,  but  this  effect  will  be  ignored,  in  order  to  avoid  con¬ 
siderable  coinplexities,  and  the  spectral  distribution  will  be  assumed  constant 
throughout  Uic  flare. 

5.2  TEMPORAL  VARIATION  OF  X-RAY  FLUX 

T'he  temporal  variation  of  the  X-ray  flu.x  during  a  flare  has  betm  estimated 
by  Poi)poff  and  Wliitten  (1961)  to  have  tlic  form  shown  on  Fig.  5.2  rising  to 
maximum  intensity  in  about  300  sec  and  then  decreasing  to  about  20  percent  of 
ils  maximum  value  in  a  .similar  interval.  Also  shown  on  Fig.  5.2  is  a  i)Iot  of 
the  fuiKdion,  Q(t),  formulated  by  Poppoff  and  Whitten  to  fit  this  variation  (this  is 
Uk.'  light  eoulinuous  line  on  Fig.  5.2). 


96 


FIG.  5.2  TEMPORAL  DEPENDENCE  OF  X-RAY  FLUX  IN  A  SOLAR  FLARE 


T’lic  statistics  j^iven  in  See.  2  suggest  tliat  tJie  decay  to  rionnal  of  the  X-ray 
flii.x  takes  rather  longer  than  is  indicated  by  Po))|)off  and  Wliitten.  Also  thcMr 
form  for  (i(t).  being  two  different  functions  in  the  intervals  0  -t’  .TOO  and  t>300, 
seems  someniiat  artificial  and  is  rather  incotivcnient  to  use  for  any  calculations 
of  tile  ioni/.ation  it  produces.  I'or  these  reasons  it  was  decided  to  find  an  ex- 
jiia-ssion  for  (^(t)  giving  a  maximum  inlen.sity  at  t  300  sec  and  a  decay  to 
(.)  Q  max.^2()  at  I  1200  sec.  A  function  of  the  form 

!i7 


(5.1) 


where  (Jo'O'.  />'  ‘'re  adjustable  con.stants  was  first  tried,  but  it  was  foiiml 

that  it  is  not  possible  to  fit  the  above  eonclitions  with  such  a  fiinetion.  The 
closest  approach  to  the  required  form  is  that  plotted  on  Fig.  .5.2,  for  o  = 

2'!)14.10  ■*  and  /I  -  :b787. 10’ *,  which  gives  C}(1200)  ~  Q  max/,5  eonsich'rably 
larger  than  is  ilesired.  The  conditions  arc  moT'C  closely  fitted  by  the  fiiiKdion 

OU)  -  (.5.2) 

with  a  =  1/150,  which  is  also  plotted  in  Fig.  5.2.  For  this  function  (^(1200) 

Q  max/25,  which  is  less  than  Q  max/20  as  required.  The  constants  and  (J,,' 
of  Eqs.  (5.1)  and  (5.2)  arc  adjusted  to  give  Q(300)  =  Q  max.  This  gives  (J,, 

Qm.ox  (0.2x10’^)  *  andQ.,'  =  Q  max  (1.22x10'’)' \ 

5.3  RESULTING  HEIGHT  AND  ZENITH  ANGLE  DEPENDENCE  OF  IONIZING 

FLUX 

The  effect  of  an  X-ray  flare  in  causing  an  SID  depends  on  the  ionization 
produced  at  any  position  in  the  ionosphere.  This  is  a  function  of  height  and 
zenith  angle  of  the  sun  at  the  time  of  the  flare  as  well  as  of  the  ionizing  flux. 
However,  the  height  and  zenith  angle  dependence  for  a  given  X-ray  flux  has 
already  been  considered  and  Figs.  3.3  through  3.9  show  the  ionization  produced 
per  erg.  cm  see  ’  of  X-ray  flux  in  terms  of  height,  zenith  angle,  and  wave¬ 
length.  By  using  these  and  the  spectral  distribution  given  in  Fig.  5.1,  the  ioniza¬ 
tion  produced  at  the  maximum  of  the  flare,  that  is,  by  a  fhtx  Q  max  lasting  for 
1  sec,  has  been  calculated,  and  the  result  plotted  on  Fig.  5.3.  By  u.sing  this  and 
a  temporal  variation  of  the  ionizing  flux  such  as  given  by  Eqs.  (5.1)  or  (5.2). 
the  ionizing  flux  during  a  moderate  flare  is  completely  specified:  that  is,  a 
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po.ssililc  I'onn  oT  Uie  behavior  ol'  the  ionizing  I'lux,  cionisislL'nl  willi  llu;  available 
experimental  observations,  is  del'inetl. 

The  slif^ht  irre^^ularities  superimposed  upon  the  curves  of  I’ig.  5.:j  aj’o 
not  of  sifvnifieanee.  In  order  to  simplify  tlic  calculations,  the  X-ray  flux,  as 
indicated  on  Fig.  5.1,  was  split  up  into  bands  centered  at  25,  16,  10,  6,  4,  2, 
and  lA.  The  irregnlarities  of  Fig.  5.3  are  the  result  of  this  subdivision. 

5.4  SOLUTION  OF  THE  BALANCE  EQUATIONS  DURING  A  FLARE 
The  equations  of  balance  for  the  ionosphere  are 


,/Y 
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where  Q(t)  is  the  ionizing  flux  during  the.  flare  and  all  other  quantities  are  as 
previou.sly  defined,  the  suffix  q  representing  quiet  conditions  before  the  flare. 
The  solution  of  these  equations  is  com|jlieatcd,  but  all  the  necessary  information 
has  been  given.  An  approximate  solution  exi^ressing  N*  as  a  step  function  with 
5-sec.  intervals  can  be  obtained  and  a  computer  program  was  set  uj)  to  do  this 
for  a  specific  example,  namely,  the  ca.se  when  the  flare  occurs  directly  over¬ 
head  at  the  equator.  The  temporal  variation  of  the  X-ray  flux  is  taken  as  given 
hy  Fq.  (5.2).  The  variation  of  N  is  tlicn  a  series  of  steps  followed  by  a  rapid 
decay  to  regain  the  attachment-detachment  balance  N  /N  A/F.  'Phis  ratio  is 
effectively  attained  in  less  than  5  sec,  so  that  before  each  addition  of  further 
N*,  the  value;  of  N  can  be  given  as  N  F  NVA+F.  This  has  also  been  e()nii)uted 
for  Hie  (;a.sc  ,\  o  and  heights  of  40  to  100  km  in  5-km  incrcmicnls,  and  Hie 


FIG.  5.4  VARIATION  IN  IONIZATION  DURING  AND  AFTER  A  SOLAR  FLARE 
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results  jiloliuti  (in  Fi^'.  5.4.  'I'lie  (iiiict  electron  densities  belure  and  alter  the 
flare  are  those  of  profile  A  in  I'igs.  4.1  and  4,2;  the  simplification  of  negleetintt 
tlic  temporal  ehan{>:e  in  solar  zenith  angle  is  adopted.  The  results  of  Fig.  4.2 
suggest  tliat  til  is  is  not  unreasonable  for  flares  occurring  witliin  ifl  hours 
(110“'  see)  of  noon.  However,  for  times  exceeding  some  lO'*  see  tlie  curves  of 
Fig.  5.4  are  unrealistic.  Interesting  features  of  Fig.  5.4  arc  that,  proportionately, 
the  maximum  added  ionization  occurs  at  a  height  of  about  85  km;  that  the  onset 
time  of  the  maximum  increases  as  height  decreases;  and  tliat,  surprisingly 
enough,  the  persistence  of  die  added  ionization  tends  to  decrease  rather  than 
increase  witli  increasing  height. 
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G.  Till-;  lUNCXSPIILUlC  DlSTTRBANCl-:  Dl  l. 


TO  A  NUCLEAR  EXPLOSION  IN  SPACE 


(i.l  INTRODUCTION 

Bolk'ii  and  Dyce  (1961)  point  out  that  most  ol'  the  energy  (60  to  SO  pereent) 
of  a  nuclear  burst  in  space  appears  in  the  form  of  prompt  X  rays„  This  raili- 
ation,  in  a  manner  very  similar  to  that  of  the  flux  of  solar  X  rays  during  an  SID. 
produces  copious  ionization  in  Ihc  lower  ionos])he.rc.  There  arc,  howevcT',  three' 
differences  between  the  disturbance  due  to  a  nuclear  event  and  that  caused  by  a 
natural  SID.  First,  and  most  important,  the  time  variations  of  the  ionizinp'  flii.x 
will  be  different;  notably,  the  time  from  onset  to  maximum  flux  will  be  almost 
Instantaneous  for  the  nuclear  burst  as  compared  with  that  for  the  natural  SID. 
Second,  the  spectral  distribution  of  the  incident  X  ray  flux  will  not  be  the-  same: 
this  will  lead  to  differences  in  the  ionization-lieight  profile  for  the  two  types  ol 
incident.  Third,  since  the  bur.st  will  not  be  coincident  with  the  sun,  the  depen¬ 
dence  upon  zenitli  angle  will  differ.  This  latter  point  is  illustrated  in  the  geo- 
mctiy  of  Fig.  G.l.  Suppose  the  burst,  N,  occurs  upon  the  line  joining  the  center 
of  the  eartli,  O,  and  that  of  the  sun,  S.  This  situation  is  that  likely  to  produce  the 
greatest  similarity  between  the  nuclear  and  the  solar  zenith-angle  d('pendenee. 
Furthermore,  for  clarification,  let  C  be  an  eijuatorial  point.  ;it  the  time'  of 
cciuinox,  experiencing  local  noon:  A  and  D  are  the  poles.  The  solar  zenith  angh'. 
V  is  effectively  equal  to  the  latitude  ",  since  the  solar  distance.  1).  is  very  miieh 
greater  than  tiu’  radius  of  tlu'  earth,  R.  'fills  is  not  true,  however,  for  the  nuc¬ 
lear  angle,  c.  .since,  for  the  most  lik(’ly  situations,  D.ll  R,  where  11  is  the  lu'ight 
of  tIu'  hurst  above  the  sniTace  of  the  earth. 

To  show  the  variation  of  zenith  angh'  with  tlu-  height  of  a  iiueb'ar  hurst,  the 
values  ol  (\  for  a  burst,  at  10^’  km  and  for  one  atio'  km  havi'  her'ii  tabulated  in 

in:! 


'ruble  C.  1  as  a  I'unetion  of  '  ,  for  the  };eotm,‘lrv  of  I'ij;-.  (i.l,  If  is  l  iiar  tliai  the 
zenith-angle  variation  from  the  solar  value  Ijecomes  apin'eeiable  bn'  cxphisinns 
fairly  close  to  the  earth. 
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G..2  SPFC'I  HAL  DIFFF.HENCES 

It  is  usual  to  assume  that  a  nuclear  burst  in  space  will  radiate  a  pulse  ol 
X  ra,\  s  whose  spi-etral  rlistribut ion  will  eoi  respond  to  that  of  a  black  body  at  a 
teniperature  ehtiraeteristie  of  the  iv\|)losion  Temperatures  of  belween  (i.  1(1  aud 
1.25  ke\'  have  been  eonsiderc-cl  by  Hollen  and  Dyei'  (IfXil).  Aei-oi'ding  to  tJie 
Wien  displacaamait  law.  \,„T  eonslant  where  A,„  is  the  wavelength  at  whiidi 
maximum  energy  is  radiat('d.  Tempeiatures  of  1.2-5  and  ('.Id  k(‘\’  eortaspond. 
res|iec:tively,  to  values  iif  -\.,,of  about  2  and  1  (( A. 

In  Fig.  (i.2  the  X-ray  intensity  incident  at  the  surfaei'  ol  the  earth  lor  the 
conditions  of  a  (|uiet  sun  and  a  moderate  Hare  a  ri'  replotted  from  tin'  infoi  ina 
tion  of  Sees.  ’■’>  and  .5.  For  eomparatixe  purposes,  the  flux  dui'  lo  a  l.GTK'T  hurst 
(1  K  T  assumed  as  energy  in  prompt  X  rays),  of  1-kev  tempiuatu I'e.  at  a  ludgld 


of  J(l"’kin,  is  alsii  shown,  lliis  ciirvo  is  ior  ihc  sub-l)Lii'si  jjoini  ill  llio  carlli,  tliiil 
is,  I'or  vcrlioiil  inc  iiloiicc'.  11  is  m  ulriit  that  thurn  are  eonsideiable  dilTereiiees 
belween  the  form  of  the  flare  eui  ve  (B)  and  that  lor  Lite  burst  (C).  As  the  si/e 
and  distaiiee  of  llie  exiilosion  are  changed,  the  over-all  magnitude  of  C  will  alter 
if  the  temperature  of  the  burst  is  modificKl,  Curve  C  will  slide  up  and  down  the 
wavelength  axis  U))on  Fig.  G.2.  Similarly,  Curve  B  will  ehange  in  size  aeeord- 
ing  to  the  magnitude  of  the  flare,  and  tliere  may  also  be  variations  in  the  spec¬ 
tral  distribution,  altliough  the  knowledge  of  these,  as  em)jhasized  in  Sec.  5,  is 
very  scanty.  However,  it  seems  certain  that  a  flare  will  always  jiriKluce  a  flux 
monotonieally  increasing  with  increasing  wavelength  over  tlie  range  of  0,.  1  to 
lOOA'  for  the  same  limits  in  wavelength  the  flax  due  to  a  burst  will  rise  to  a 
maximum  (probably  at  a  wavelength  between  1  and  loA)  and  then  decrease  us  tlu 
wavelength  further  increases.  Ignoring  the  temporal  variation,  although  it  may 
be  conceivable  for  the  X-ray  flux  from  nuclear  and  solar  disturbances  to  be 


FIG.  6.2  X-RAY  FLUX  FOR  THREE  DIFFERENT  TYPES  OF  SOURCE 
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similar  over  a  wide  range  of  wavelengths  (for  example,  by  an  effective  super¬ 
position  of  the  lower  parts  of  Curves  B  and  C),  it  seems  impossible  for  this 
similarity  to  extend  over  the  whole  wavelength  band  from  0.1  to  lOOA.  These 
differences  in  spectral  distribution  inevitably  have  ionospheric  consequences 
which  should  be  detectable. 

As  an  example  of  the  ionization  differences  that  can  result  from  a  spectral 
variation,  tiie  electron  densities  produced  by  the  nuclear  flux  (Curve  C  of  Fig. 

6.2)  and  during  1  sec  by  the  flare-associated  flux  (Curve  B  of  Fig.  6.2),  assum¬ 
ing  no  losses  by  recombination  or  attachment,  are  plotted  on  Fig.  6.3.  Both 
curves  give  maximum  electron  densities  of  the  order  of  10^  for  heights  below 
100  km,  but  the  height  of  the  maximum  for  the  specific  nuclear  example  is  much 
less  than  that  for  the  flare-produced  ionization.  This  could  have  widely  differ¬ 
ent  consequences  upon  radio  wave  propagation.  For  instance,  VLF  propagation 
is  controlled  during  daytime  by  an  electron  density  of  about  10^  cm"®  at  a 
height  of  some  70  km.  The  added  ionization  in  a  second  at  70  km  due  to  the 
flare  is  small  compared  with  10^;  thus  the  effects  on  VLF  propagation  would  be 
minor.  However,  for  the  nuclear  instance  the  enhancement  in  ionization  at  70  km 
is  comparable  with  that  already  present;  in  consequence,  the  modifications  in 
VLF  propagation  would  be  very  considerable. 

The  difference  in  spectral  distribution  of  the  flux  due  to  a  solar  flare  and 
that  due  to  a  nuclear  burst  will  also  result  in  considerable  differences  in  zenith- 
angle  dependence  of  the  resulting  ionization.  A  height  of  10^  km.  as  used  by 
Bollen  and  Dyce  (1961)  for  the  nuclear  burst,  has  been  considered,  for  a 
1.G7KT  nuclear  burst  at  1-kev  temperature,  which  gives  a  1-KT  yield  of  X  rays 
having  the  spectral  distribution  shown  on  Curve  C  of  Fig.  6.2.  Using  this  curve 
and  interpolating  from  the  curves  of  Figs.  3.3  through  3.8.  the  ionization 
produced  at  the  angles,  a,  given  on  Table  6.1,  corresponding  to  e  =  0°,  20°,  40°, 
60°,  and  80°  has  been  calculated  for  heights  from  40  to  100  km  and  the  results 
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plotted  on  Fig.  6.4.  The  irregularities  of  Fig.  6.4  are  probably  introduced  Iqr 
approximating  the  continuous  X-ray  spectrum  by  a  definite  flux  at  each  of 
several  discrete  wavelengths  (the  same  is  true  of  the  calculations  of  the 
ionization  due  to  a  solar  flare  shown  on  Figs.  5.3  and  6.3).  A  comparison  of 
the  curves  of  Fig.  6.4  with  those  of  Fig.  5.3  shows  considerable  differences  in 
over-all  behavior,  in  both  the  zenith-angle  and  height  dependence  of  the 
ionization. 

It  may  be  concluded  that,  even  if  the  X-ray  flux  during  an  unusually  abrupt 
solar  flare  occurred  in  a  burst  of  magnitude  and  duration,  as  observed  by 
detectors  at  the  earth,  comparable  with  that  associated  with  a  nuclear  burst  in 
space,  then  differences  in  the  spectral  distribution  of  the  X  radiation  would  pro¬ 
duce  corresponding  ionospheric  variations  which  could  be  distinguished  by 
radio  methods. 

6.3  THE  DECAY  OF  IONIZATION  FOLLOWING  A  NUCLEAR  BURST  IN  SPACE 

If  the  pre-existing  ionization  is  neglected,  and  the  deposition  of  ionization  by 
the  nuclear  burst  taken  as  instantaneous,  with  zero  production  before  and  after 
the  explosion,  then  the  ionospheric  balance  equations  with  B  =  E  are  of  a  very 
simple  form: 


—  =  _ /?;V,V  +  FV 

dt 


^  =  A^  -  FN-  -  RN'N* 
dt 
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with  the  solutions 
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where  N,,  is  the  positive- ion  or  electron  density  produced  at  time  t  =  O. 

The  variation  of  N  with  time  is  plotted  in  Figs.  6.5  through  6.10  for  values 
of  Ng  of  10*,  10*,  10'*,  and  10* ,  and  heights  at  10-km  intervals  between  60  and 
100  km.  Both  day  and  night  circumstances  are  shown.  X  is  taken  as  0°  during 
the  day.  It  is  evident  that  the  decay  of  the  electron  density  often  has  three  main 
phases.  Initially,  the  loss  of  electrons  is  primarily  by  attachment,  until  a  quasi¬ 
balance  of  attachment  and  detachment  processes  has  been  established,  with  a 
corresponding  partition,  almost  independent  of  N^,  of  negatively  charged  part¬ 
icles  between  electrons  and  negative  ions.  Subsequently,  there  is  a  plateau 
during  which  the  electron  density  remains  approximately  constant;  the  duration 
of  the  plateau  depends  upon  the  height  and  upon  N^.  Finally,  the  decay  phase  of 
the  negative  and  positive  carriers  by  recombination  takes  over.  There  is  then 
a  strong  tendency  for  the  curves  for  different  values  of  to  approach  each 
other;  in  other  words,  the  remaining  electron  density  after  a  very  long  time  is 
independent  of  N^.  It  is  noteworthy  that  the  duration  of  the  excess  ionization 
does  not  increase  entirely  systematically  as  the  height  increases;  this  is  exem¬ 
plified  on  Fig.  6.10  for  the  night  curves. 
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FIG.  6.8  DECAY  OF  NUCLEAR  IONIZATION  -  HEIGHT  90  km 
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FIG.  6.9  DECAY  OF  NUCLEAR  IONIZATION  -  HEIGHT  100  km 
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FIG.  6,10  DECAY  OF  NUCLEAR  IONIZATION  -  Np  -  10^ 
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The  treatment  is  altered  somewhat  if,  with  all  other  conditions  unchanged, 
tlie  effects  of  a  constant  rate  of  production  of  ionization  before  and  after  the  ex¬ 
plosion  and  the  resulting  equilibrium  value  of  the  ionization  are  included.  The 
appropriate  equations  are 


rf.Y 

=  o  _  ,gv  + /'Y- - /n',v* 

(It 


^  ^  ,1/V  - 
(It 


=  Q- 
(It 


with  the  initial  conditions  that  N  =  Nq  +  N^,  N*  =  N'  =  N  "  at  time 

t  =  o.  Suffix  q  gives  the  equilibrium  values  for  the  quiet-sun  conditions  which 
were  considered  in  Sec.  3.  By  writing  No  =  KN^  and  considering,  for  example, 
heights  of  60-70  km,  it  is  found  that  for  K  <100  it  is  possible  to  solve  these  equa¬ 
tions  approximately  by  considering  two  separate  stages  in  the  return  to  equili¬ 
brium  conditions,  as  suggested  by  Swift  (1961).  In  the  first  stage,  tlie  effects  of 
recombination  and  the  steady  ionization  are  neglected  so  that  remains  sensi¬ 
bly  constant  while  there  is  a  rapidly  achieved  balance  of  attachment  and  detach¬ 
ment  resulting  in  the  ration  N*/N  =  A/F,  which  is  then  maintained.  The  second 
stage  consists  of  a  slow  decay  of  electrification  to  the  equilibrium  values,  by  re¬ 
combination.  The  solutions  are  then,  for  the  first  stage 
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and  /V*  =  ,V*-A'; 

and  for  the  second  stage 
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For  the  above  treatment  to  be  valid,  it  is  necessary  tliat  N^docs  not 
change  appreciably  during  the  short  time  during  which  the  ratio  N"/N  is  becom¬ 
ing  insignificantly  different  from  A/F.  Two  conditions  are  involved.  The  first 
is  that 
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which  gives 
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The  second  condition  is  satisfied  when 
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which  is  obeyed  for  t  such  that 


%(A  +  F)2 


These  together  require  K  to  be  such  that 


AK  \  3  (/I  +  F)^ 

A  +  (a:  +  l)F  I  2,/HQ  {k^F^  +  (/I  +  F)^) 


With  K<100  and  values  of  the  parameters  A,  F,  B,  and  Q  appropriate  for  heights 
of  60  to  70  km,  the  above  inequality  is  certainly  satisfied. 

The  solutions  for  K  =  1  and  K  =  10  at  heights  of  60,  65,  and  70  km  have 
been  plotted  in  Figs.  6.11  and  6.12  using  the  values  of  the  constants  for  zenith 
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angle  and  magnetic  dip  angle  both  zero,  that  is,  neglecting  any  variation  in  Q 
and  F  due  to  changing  zenith  angle  during  the  time  of  the  return  to  equilibrium. 
For  comparison,  the  curve  obtained  neglecting  Q  for  a  similar  initial  value  of 
total  ionization  has  also  been  drawn  (Fig.  6.11)  corresponding  to  the  case  K  =  10 
and  a  height  of  70  km.  It  is  seen  to  be  almost  identical  in  the  initial  stages,  but 
differs  greatly  as  t  increases  so  that  the  effect  of  the  steady  Q  begins  to  show. 
The  differences  are  in  the  later  recombination  stage  rather  than  in  the  earlier 
stage  of  partition  of  the  negative  carriers  between  electrons  and  negative  ions. 
Further  calculations  showing  the  importance  of  considering  pre-existing  ioni¬ 
zation  are  included  in  Appendix  I,  in  which  the  relevance  of  the  present  work  to 
two  recent  reports  by  Westfall  and  by  LeLevier  is  considered. 


FIG.  6.11  DECAY  OF  NUCLEAR  IONIZATION  WHEN  PRE-EXISTING  IONIZATION 
IS  CONSIDERED  (K  -  10) 
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FIG.  6.12  DECAY  OF  NUCLEAR  IONIZATION  WHEN  PRE-EXISTING  IONIZATION 
IS  CONSIDERED  (K  *  11 
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7.  SUMMARY  AND  CONCLUSIONS 


This  report  is  primarily  designed  to  provide  a  survey  of  the  phenomenology 
of  SEAs  and  a  background  of  ionospheric  tlieory  for  use  in  further  work.  It  is 
not  intendwl  to  be  a  self-contained  and  complete  entity.  However,  it  is  still 
desirable  to  summarize  the  major  points  that  have  emerged.  These  arc: 

(1)  There  is  an  association  but  no  rigid  interconnection  between  the  various 
SID  phenomena  accompanying  a  solar  flare. 

(2)  SEAs  have  an  average  rise  time  of  about  10  min.  Rise  times  of  less 
than  2  min  are  relatively  rare,  but  occur  upon  the  average  at  a  rate 
of  perhaps  one  every  two  months. 

(3)  SEAs  due,  respectively,  to  a  solar  disturbance  and  to  a  local  thunder¬ 
storm  have  considerable  similarity.  The  best  way  of  differentiation 
is  to  use  the  fact  that  an  Sf2A  produced  by  thundery  activity  gives  an 
enhancement  at  all  frequencies;  the  SEA  accompanying  an  SID  has 
an  enhancement  at  frequencies  exceeding  20  kc  and  a  decrease  in 
noise  for  frequencies  less  than  10  kc. 

(4)  There  are  considerable  differences  in  the  temporal  behavior  during 

an  SEA  at  frequencies  of  10  and  27  kc.  It  seems  probable  that  the  10-kc 
noise  is  influenced  by  ionization  at  a  lower  level  than  is  the  case  for 
a  frequency  of  27  kc.  This  implies  that  the  significant  ionizing  radi¬ 
ations  and  the  speed  of  ionospheric  recovery  are  different  for  the  two 
frequencies. 

(5)  The  interesting  behavior  as  a  function  of  frequency  during  an  SEA 
suggests  that  simultaneous  observations  at  many  frc-cjuencies  in  the 
VLF  band  would  be  valuable. 

((!)  A  survey  of  information  on  the  rate  coefficients  for  the  lower  ionos|iherc 
suggests  that  the  values  of  tlie  att:ichmcnt  coefficient  arc-  wc-ll  established. 
As  regards  detaelimont,  photodetachment  is  dominant  by  day  at  heights 
exccc'ding  (it)  km;  below  tlii.s  level,  collistonal  dehiehmcnt  seems  signif¬ 
icant.  At  night  collisional  detachment  dominates  below  70  km  while 
associative  detachment  is  jirobably  important  at  greater  heights.  The- 
elc'ctron-ion  and  ion- ion  recombination  coc-ffieients  are  of  comiiarable 
magnitude  and  do  not  vary  mucdi  with  height.  Consc-qiicntly,  tlie  eonsid- 
c  rable  simplification  in  ionospheric  tliooi'y  of  regarding  the  two  recom¬ 
bination  eoeffic-.ient.s  as  equal  and  constant  seems  valid. 


(7)  TIic  iiiajoi'  ionixin*;  af^iiicics  in  llio  lower  itjnosplierc  are  .sohii'  X  rays, 
solar  l.yinan  a  ,  and  eosmu:  rays.  The  latter  is  die  only  signilieanl  in- 
lliieiicc  at  ni;_’;hl,  ami  even  liy  tlay  is  cloniinunt  below  (iO  Uni.  In  daytime 
and  tor  hei^;hts  Iietween  GO  and  80  km,  tlie  three  ageneies  arc  ol'  com¬ 
parable  strengih;  above  80  km,  X  rays  are  most  inpiortant.  These 
conclusions  are  for  a  period  near  solar  maximum;  lliere  are  indications 
that  the  significance  of  solar  X  rays  relative  to  solar  Lyman  a  may 
decrease  with  sunspot  number. 

(8)  The  balance  equations  for  the  lower  ionosphere  can  be  solved  easily, 
assuming  no  vertical  transport  of  ionization  and  cciuating  the  recombi¬ 
nation  coefficients.  By  day  tlie  equilibrium  condition  corresponding  to 
a  very  slow  rate  of  change  of  ion  density  is  closely  appioached;  at 
night,  however,  decay  of  the  ionization  from  an  equilibrium  daytime 
value  must  be  considered.  There  arc  considerable  global  asymmetries 
in  electron  density  contours.  Generally  speaking,  however,  two  electron 
layers  can  be  identified  at  night;  one  centered  at  about  55  km  is  due  to 
cosmic  rays,  while  alxjvc  some  80  km  ionization  produced  during  the 
day  persists  throughout  the  night.  With  the  onset  of  daytime  conditions, 
the  ionization  above  (>5  km  increases  rapidly. 

(9)  A  parameter,  cUy,  which  may  be  considered  as  an  effective  ionospheric 
conductivity,  defines  the  upper  boundary  of  the  earth-ionosphere  wave¬ 
guide  system  involved  in  VLF  propagation.  Conversion  of  electron 
density  contours  into  contours  of  shows  that  these  are  not  parallel 
to  the  surface  of  the  earth,  particularly  l)y  day.  For  meridional  propa¬ 
gation  the  contours  of  co,.  are  vciy  different  according  to  whether  the 
term  in  co,.  involving  the  earth's  magnetic  field  is  included  or  neglected. 

(10)  During  a  moderate  solar  flare  the  flux  of  X  rays  affecting  the  lower 
ionosphere  increases  by  a  factor  of  pcrhajis  twenty.  The  flux  rises  to 
its  maximum  in  a  time  of  some  5  minutes;  the  subsequent  decay  occupie 
about  15  minutes.  The  X-ray  flux  produces  its  greatest  proportionate 
increase  in  ionization  (by  a  factor  of  about  four)  at  a  height  near  80  km 
The  time  to  reach  maximum  electron  density  varies  between  300  and 
1000  sec.  according  to  height;  the  greatest  lag  is  at  the  lower  heights. 
Suiqirisingly  enough,  the  rate  of  decay  of  the  added  ionization  tends  to 
increase  with  increasing  height  . 

(11)  A  nucleai'  burst  in  space  emits  a  flux  of  X  rays  iliat  differs  in  thiaa' 
significant  rcs|>ecls  from  the  flux  associated  with  a  solar  flai’e. 

First, — and  most  im|Jortanl — the  onset  is  immediate.  Second,  the 
spectrum  of  th<'  nuclear  flux  is  conccnti'ated  more  towards  energetic 

X  rays  than  is  the  solar  flux.  Third,  the  zenitti-angte  dependence's  arc 
different,  't  he  more  encrge'tic.  spectrum  implies  the  production  of 
ioniz.ation  at  low  leve'ls  in  tiie  ionosphere.  Coml)ined  witli  the  sudden 
onset,  this  means  that  tlie  ionization  changes  at  a  height  of  70  km,  a 
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controlling  level  in  daytime  VLF  propag'ation,  are  far  more  I’apicl  during 
a  nuclear  burst  than  during  a  solar  flare. 

(12)  After  the  sudden  deposition  of  ionization  in  the  lower  ionosphere  by  a 
nuclear  burst  in  space,  three  decay  phases  can  often  be  identified:  the 
establishment  of  an  initial  partition  of  negative  carriers  between  elec¬ 
trons  and  negative  ions;  a  plateau  stage  during  which  there  is  little 
change  in  electron  density;  and  a  final  phase  in  which  recombination 
processes  control  the  decay.  The  pre-existing  ionization  has  consider¬ 
able  influence  upon  the  decay,  especially  during  the  recombination  phase. 

(13)  If  a  specific  electron  density  is  suddenly  deposited,  as  by  a  nuclear 
burst,  at  different  levels  in  the  ionosphere,  there  are  apparent  anoma¬ 
lies  in  the  decay  of  the  ionization.  When  account  is  taken  of  the  pre¬ 
existing  ionization,  the  rate  of  decay  does  not  decrease  systematically 
as  height  increases,  as  might  be  anticipated.  The  exact  specification 
of  this  anomalous  behavior  depends  critically  upon  the  values  assumed 
for  the  rate  coefficients  and  for  the  quiet-state  ion  production.  An 
example,  for  heights  of  70,  75,  and  80  km,  shows  that  in  the  attachment- 
detachment  phase  electrons  are  lost  more  rapidly  at  the  lower  heights. 
The  contrary  is  the  case  for  the  recombination  phase,  so  that  ultimately 
the  greatest  electron  density  exists  at  70  km. 

(14)  It  is  shown  that  the  departures  from  parallelism  of  the  earth  and 
ionosphere  boundaries  for  the  waveguide  theory  of  VLF  propagation 
have  significant  Implications.  The  true  attenuation  coefficient  can 
differ  from  an  average  value,  deduced  for  an  ionospheric  boundary 

of  constant  height,  by  10  percent  or  more  according  to  the  time  of  day. 
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REPORTS  BY  WESTFALL  AND  LELEVIER 


I.l  GENERAL 

Since  the  completion  of  much  of  the  work  of  this  report,  two  important 
relevant  reports  have  become  available,  and  it  is  worthwhile  comparing  their 
results  with  those  already  stated.  The  first  of  these,  by  Westfall  (1962)  deals 
with  the  variation  of  the  electron  density  in  the  ionosphere  caused  by  an  SID,  in 
particular  as  observed  in  the  sudden  phase  anomaly  (SPA)  occurring  in  VLF 
propagation.  From  observations  of  this  phenomenon  it  is  possible  to  calculate 
the  changes  in  the  effective  height  of  reflection  for  a  given  frequency,  which  can 
then  be  related  to  electron  density  at  that  height.  Two  points  arise  that  are  of 
interest  here.  The  first  is  that,  as  for  the  SCNA  and  SEA  statistics  considered 
in  Sec.  2,  the  SPA  observations  did  not  indicate  that  there  is  any  close  correlation 
between  the  visible  Ha-flare  and  the  ionizing  flux,  either  in  temporal  behavior 
or  magnitude.  A  lack  of  correlation  is  also  noted  in  comparing  the  ionizing  flux 
to  2800- Me  solar  noise  bursts,  although  all  three  phenomena  are  generally  re¬ 
lated  and  are  products  of  a  disturbed  state  of  the  sun.  Second,  Westfall's  de¬ 
ductions  of  spectral  behavior  during  a  flare  are  worth  noting.  The  relation 
between  X-ray  flux  and  electron  density  produced  is  considered  to  be  that  given 
by  Pierce  and  Swift  (1960).  The  rate  coefficients  used  in  that  report  are  con¬ 
siderably  different  from  those  given  in  Sec.  3  of  this  report,  but  the  same  gen¬ 
eral  behavior,  a  nonlinear  decrease  in  electron  production  witli  increasing  solar 
zenith  angle  and  decreasing  ionospheric  height,  is  found  witli  the  modified  para¬ 
meters  as  was  obtained  with  the  Pierce  and  Swift  (1960)  values.  The  solar 
zenith-angle  dependence  of  the  electron  density  caused  by  an  SID  is  largely  re- 
Intixl  to  the  spi;ctral  content  of  a  flare.  Because  the  relationship  between  solar 
zenith-angle  difference  and  the  difference  in  change  of  height  of  reflection 


129 


between  two  ])aths  during  a  flare  varies  for  different  flares,  Westfall  concludes 
that  the  spectral  content  of  tlic  X-ray  flux  is  variable.  The  treatment  of  Sec.  5 
assumes,  for  simplicity,  a  specific  spectral  distribution  although  stating  that 
this  is  not  necessarily  maintained,  even  during  tlie  course  of  a  single  flare. 
Variations  in  spectral  distribution  between  individual  flares  are  certainly  to  be 
anticipated  and  would  account  for  Westfall's  observations.  There  are,  however, 
other  factors  involved  which  could  also  be  considered  in  explaining  the  varia¬ 
tions  noted.  For  example,  the  actual  zenith  angle  and  the  latitude  of  the  path  will 
cause  differences  in  the  original  pre-flare  reflection  height  over  that  path,  and 
the  change  in  height  caused  by  a  given  X-ray  flux  is  height  dependent. 

A  second  and  more  closely  related  report  is  that  of  R.  E.  LeLevier  (1962) 
which  covers  essentially  the  material  of  Sec.  6.3.  A  summary  of  the  rate  co¬ 
efficients  involved  in  determining  ionospheric  balance  conditions  is  made,  but 
it  is  not  as  detailed  as  that  of  Sec.  3.  A  brief  comparison  of  the  values  obtained 
is  made  here. 

1.2  ATTACHMENT  (A) 

LeLevier  considers  A  to  be  solely  due  to  a  three-body  attachment  process 
whereas  in  Sec.  3  it  is  found  that,  above  about  80  km,  radiative  attachment  to 
Og  and  also  attachment  to  atomic  oxygen  become  significant.  Figure  1.1  gives 
an  approximate  comparison  of  the  two  sets  of  values. 

1.3  ELECTRON-ION  RECOMBINATION  (B) 

LeLevier  uses  the  value  of  4  x  10"’,  corresponding  to  recombination  with 
02*  ,  which  is  approximately  the  upper  bound  given  in  Sec.  3.  This  is  an  order  of 
ten  greater  than  that  used  Oiroughout  this  report  (B  =  4  x  10'®  ),  a  difference 
which  considerably  modifies  the  results  obtained  when  considering  the  decay  of 
ionization.  The  reasons  for  the  value  employed  in  this  report  are  given  in  Sec.  3. 
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FIG.  1.1  COMPARISON  OF  VALUES  OF  PIERCE  AND  LE  LEVIER  FOR  THE  ATTACHMENT 
COEFFICIENT 
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1.4  ION-ION  RECOMBINATION  {E) 


In  most  of  tlie  calculations  of  this  report  E  is  taken  as  equal  to  B,  and  is 
therefore  given  by  4  x  10"“.  LeLevier  carries  out  his  analysis  twice,  in  one 
case  equating  E  to  B,  and  in  the  otlier  using  E  =  O.  However  the  difference  in  B 
already  noted  results  in  a  different  E  even  when  the  same  assumption  (i.e.,  E  = 

B)  is  made. 

1.5  DETACHMENT  (F) 

In  this  report  F  =  C  +  D  where  C  is  the  term  due  to  photodetachment  and  D 
is  the  sum  of  collisional  and  associative  detachment.  The  figures  given  show  D 
to  be  significant  below  about  70  km,  whereas  LeLevier  considers  it  is  negligible 
and  includes  only  photodetachment  from  the  O2'  ion.  This  differs  even  from,  the 
C  term  of  this  report,  because  above  80  km  the  contribution  of  detachment  from 
the  O' ion  is  appreciable.  LeLevier  makes  no  mention  of  zenith-angle  dependence, 
implicitly  assuming  a  zenith  angle  of  zero  in  his  value  of  the  photodetachment 
coefficient. 

The  values  of  the  ambient  electron  density  (N^)  have  been  shown  in  Sec.  3 
to  vary  with  latitude  and  zenith  angle  as  well  as  with  height.  In  Fig.  1.2  LeLevier 's 
assumed  electron  density  profile  is  compared  with  that  obtained  for  quiet-sun 
conditions  and  latitude  and  zenith  angle  both  zero. 

LeLevier  obtains  an  exact  solution  for  the  decay  of  ionization  after  a  sudden 
deposition  of  ions  which  gives  the  same  result  as  the  approximate  solution  of 
Sec.  6.3.  This  is  true  during  both  stages  of  the  decay,  in  all  the  cases  con¬ 
sidered,  using  either  set  of  parameters,  since  the  most  important  condition, 

VBQ<<  A  +  F,  holds.  Thus  any  differences  in  results  obtained  arc  due  to  the 
different  choice  of  values  for  the  parameters.  Choosing  the  height  67  km,  at 
which  the  ambient  electron  densities  of  Fig.  1.2  are  equal,  tlie  decay  of  ioniza¬ 
tion  after  a  deposition  of  10^  elcctrons/cm^  has  been  calculated  using  both  sets 
of  parameters.  The  excess  ionization  (N-N^)  is  plotted  on  logarithmic  scale  in 
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10 


FIG.  1.2  CO.MPARISON  OF  QUIET  SUN  ELECTRON  DENSITY  PROFILES  AS  GIVEN  BY 
PIERCE  AND  BY  LE  LEVIER 
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Fiy,iirc  1.3,  Figure  1.4  shows  the  total  electron  densities  for  eomparison  with 
earlier  figures  in  tliis  report. 

One  of  the  most  interesting  features  of  LeLevier's  work  is  his  confirmation 
of  an  effect  originally  noted  by  BocKer  and  Crain  (1962).  This  work  showed  that 
if  a  certain  density  of  electrons  is  suddenly  added  at  different  levels  in  the 
lower  ionosphere,  then  the  rate  of  decay  of  the  additional  ionization  docs  not 
decrease  steadily  as  height  increases,  as  might  be  expected.  There  are  regions 
where  the  persistence  of  the  ionization  is  greater  than  the  levels  above  or  below. 
Booker  and  Crain  found  that  these  apparently  anomalous  persistence  effects  are 
especially  marked  when  account  is  taken — as  it  obviously  should  be — of  any  pre¬ 
existing  steady  ionization.  Figure  6,10  (night)  shows,  however,  that  these 
anomalies  can  be  present  even  when  the  pre-existing  ionization  is  ignored. 

For  daytime  conditions  and  with  the  parameter  B  =  E  =  4  x  10'’  LeLcvier 
gives  curves  showing  that  the  decay  of  an  added  10^  electron/ cm^  is  less  rapid 
at  75  km  than  it  is  at  70  or  80  km,  A  similar  analysis  has  been  performed  with 
the  parameters  of  the  present  report  (notably  B  =  E'  =  4  x  10"®  )  and  the  results 
are  shown  in  Fig.  1.5,  those  of  LeLevier  also  being  given  for  comparison.  Dur¬ 
ing  the  phase  in  which  the  attachment-detachment  balance  is  being  established, 
that  is  up  to  between  perhaps  10  and  100  sec,  both  sets  of  results  show  that  the 
proportion  of  the  ionization  persisting  increases  with  increasing  height.  How¬ 
ever,  for  times  between  10^  and  10“*  sec  when  the  ion  decay  phase  has  taken 
over,  LeLevier's  parameters,  as  already  mentioned,  lead  to  a  greater  electron 
density  at  75  km  than  at  70  or  80  km.  When  the  Pierce  parameters  are  em¬ 
ployed,  the  even  more  surprising  result  is  obtained  that  the  persisting  electron 
density  in  the  ion  decay  phase  becomes  greater  at  70  km  than  it  is  at  either  75 
km  or  80  km. 

It  is  evident  tliat  tliese  apparently  anomalous  persistence  effects  are  quite 
common.  Also,  it  is  difficult  to  establish  their  exact  magnitude  and  significance 
without  precise  knowledge  of  such  quantities  as  the  ambient  electron  density, 


the  quiet-slate  ion  pi'oduction,  and--most  iniportant--tlie  rate  coefficients. 
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FIG.  i.5  DECAY  OF  DFPOr,irF  D  11  ECTRONS  -  COMPARISON  OF  LE  LEVIER 
AND  PIERCE  PARAMETERS 
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Appendix  II 


THE  EFFECT  OF  VARIABLE  IONOSPHERIC  HEIGHT  ON  VLF  PROPAGATION 


The  waveguide  mode  tlioory  of  VLF  propagation  has  been  largely  formulated 
in  terms  of  an  ionosphere  parallel  to  the  surfaee  of  the  earth.  Comparatively 
little  eonsideration  has  been  given  to  analysis  of  the  influenee  of  variable  iono¬ 
spheric  height,  although  Wait  (1961)  has  examined  the  influence  of  a  local  pertur¬ 
bation.  Wait  (1961’)  has  also  treated  the  case  of  transition  from  day  to  night  prop¬ 
agation,  by  conceiving  the  problem  to  be  of  propagation  in  two  parallel-plate 
waveguide  regions  connected  by  a  tapered  section. 

It  is  evident  from  the  results  deduced  in  Secs.  4,  5,  and  6  of  this  report  that 
contours  of  equal  electron  density  in  tlie  lower  ionosphere  are  not  parallel  to  the 
earth's  surface.  This  is  true  both  for  quiet  conditions  and  also  svhen  the  ionosi)here 
is  disturbed  either  by  natural  or  nuclear  events.  The  deviations  from  parallelism 
are  comparatively  gradual  except  near  a  sunrise  zone.  Nevertheless,  svhen  prop¬ 
agation  over  great  distances  is  involved,  as  is  commonly  the  case  at  VLF,  the 
deviations  certainly  have  appreciable  influence. 

A  specific  value  of  the  parameter  Wr  ,  which  is  directly  related  to  the  ionos¬ 
pheric  conductivity,  is  convenient  in  defining  the  upper  boundary  of  the  waveguide. 
Suppose  u)r  is  chosen  to  be  10^  cps.  Tlien  the  earth- ionosphere  guide  is  taken  to 
be  sharply  bounded  at  the  height  for  which  ~  10^-.  above  this  height  the  ionos¬ 
pheric  conductivity  is  assumed  to  be  uniform,  while  below  the  level  of  lOj.  --  10®, 
cj,.  is  taken  to  be  zero.  In  order  to  illustrate  the  variation  in  height  of  tlie  contour 

ig.  II. 1  has  been  constmicted.  This  shows  the  level  uj,  “  10®  as  a  function 
of  zenith  angle  (relative  longitude)  for  propagation  around  the  equator  at  an  equinox 
and  under  quiet  ionospheric  conditions.  'I'he  re.sults  used  in  preparing  Fig.  ILl 
were  derived  in  Sec  .  4.  It  is  notewortliy  that  tlie  c'urvc  for  height  versus  zenith 
angle  is  asymmetric,  with  the  greatest  slope  at  sunrise. 
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FIG.  il.l  VARIATION  OF  EFFECTIVE  IONOSPHERIC  HEIGHT  WITH  TIME  OF  DAY 
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Wait  Irc-atod  a  iocal  perturbation  by  postulating  a  cosine  iunction  lor 

the  variation  ol  ionosplioric  height  witli  liori/.oiital  distance.  The  same  type  of 
approa(;h  can  legitimately  be  applied  over  iniieli  greater  distances  to  obtain 
order  of  magnitude  estimates  for  the  iiiflueiiee  of  tlie  ionospheric  lieight  varia- 
.  'lily.  Some  calculations  have  been  made  for  the  height  variation  rejircsented  by 


1,(1))  ^ 


(II. 1) 


where  e  is  the  zenitlt  angle  (relative  longitude)  in  degrees.  This  eurve  is  plotted 
on  Fig.  II.  1  for  the  values  of  tlie  constants  a  =  73  and  b  =  6.5,  h  (d),  the  height, 
being  in  kilometers.  In  the  mode  theory  of  VLF  propagation,  the  major  contri¬ 
bution  to  the  deerease  in  amplitude  with  distanec  in  a  mode  arises  from  a  factor 
exp  (-ax),  where  x  is  die  distance  and  a  tlie  attenuation  coefficient.  For  the 
least  attenuated  mode  (n  =  1),  the  attenuation  factor,  a,  is  given  approximately 
(Wait  1957)  by 


a 


A„\ 


(n.2) 


where  cJis  tlie  angular  frequency  being  propagated,  A  the  corresponding  wave¬ 
length,  and  h  tlie  ionospheric  height. 


■ 


Wlien  Uie  ionospheric  height  is  variable,  Efj.  (II. 2)  must  be  modified  by  die 
inclusion  of  an  integral.  The  attenuation  factor  between  x  and  is  then  given  by 


(■’'2  dx 
J,  Ui(x)) 


(11.3) 


If  the  height  variation  of  Eq.  (II. 1)  is  adopted,  the  integral  can  be  evaluated  in 
closed  form  (x  -  oR,  where  R  is  the  radius  of  the  earth).  The  results  obtained 
at  a  frequency  of  10  ke  are  given  in  Table  II. 1. 


Table  II.  1 

INFLUENCE  OF  VARIABLE  IONOSPHERIC  HEIGHT 
ON  ATTENUATION  COEFFICIENT 


0^  (degrees) 

02  (degrees) 

h,(km) 

(km) 

a  (db/ 1000  km) 

0 

40 

66.5 

68.4 

2.34 

40 

60 

68.4 

70.5 

2.24 

60 

80 

70.5 

73.1 

l.!)l 

80 

100 

73.1 

75.7 

1.75 

100 

120 

75.7 

77.6 

1.58 

120 

140 

77.6 

70.0 

1.44 

For  a  height  of  70  km,  which  is  that  eommonly  adopted  in  daytime  VLF  propaga¬ 
tion,  the  attenuation  coefficient  as  derived  from  F.q.  (II. 2)  is  2.08  db/1000  Ion. 

A  comparison  of  this  value  with  die  figi'res  in  Table  il.l  shows  that  the  attenua¬ 
tion  coefficient  ean  differ  according  to  the  time  of  day  by  some  .+10  percent  from 
the  average  daytime  value.  The  same  re.sult  ean  be  derived  very  simply  from 
Eep  (II. 2).  Since  the  lO’  level  shown  in  Fig.  11,1  is  about  68.5  km  at  noon 

and  nc'ar  71  km  at  16.00  (  o  "  Mf).  it  follows  from  Eq.  (TI.2)  that  the  ratio  of 
attc^nuation  c-oefficients  at  tliese  times  i.s  (7  1/68.5)”  ~  l.I. 


Ml 


Obviously  tlie  asymmetric  curve  oi  Fig'.  11.1  for  a  level  of  constant  can 
be  better  approximated  by  expressions  other  than  tliat  of  Eq.  (II.  1).  However, 
the  above  analysis  is  perhaps  sufficient  to  indicate  the  potential  errors  implicit 
ill  assuming  the  attenuation  coefficient  to  be  constant  throughout  the  day  (or 
night,  for  that  matter).  Certainly  with  a  long  path,  for  which  tlie  solar  zenith 
angle  is  significantly  different  from  the  average  value  over  the  day,  the  use  of 
a  mean  daytime  attenuation  coefficient  can  lead  to  appreciable  inaccuracies. 
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